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Abstract 

This paper is devoted to the two-loop QED radiative corrections to the decay ir° — > e + e~. 
We compute the virtual corrections without using any approximation and we take into ac- 
count all the relevant graphs with the inclusion of those omitted in the previous approximative 
calculations. The bremsstrahlung is then treated within the soft photon approximation. We 
concentrate on the technical aspects of the calculation and discuss in detail the UV renor- 
malization and the treatment of IR divergences within the dimensional regularization. As a 
result we obtain the 0(a 3 p 2 ) contribution in closed analytic form. We compare the exact two- 
loop results with existing approximative calculations of QED corrections and find significant 
disagreement in the kinematical region relevant for the KTeV experiment. 
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1 Introduction 

The rare decay of the neutral pion into the electron-positron pair provides an interesting tool to 
test the nonperturbative low-energy dynamics of the Standard Model (SM). While the possible 
contributions of the weak sector of the SM are tiny and can be safely neglected, the leading order 
QED contribution is described by two virtual photon exchange diagram and is therefore tightly 
connected to the doubly off-shell pion transition form factor i ? 7r o 7 * 7 * for the subprocess 7r° — > 7*7*. 
Better understanding of this form factor which is not known from the first principles is important 
e.g. for the determination of the light-by-light hadronic contribution to the muon anomalous 
magnetic moment g — 2. On the other hand the rareness of the decay which is suppressed with 
respect to the 7r° — > 77 decay by a factor of 2(am/M n o) 2 within the SM (here m is the electron 
mass which enters here as a consequence of the approximate helicity conservation) makes it also a 
promising process possibly sensitive to the physics beyond the SM. 

The systematical theoretical treatment of the process dates back to 1959 when the first pre- 
diction of the decay rate [1] was published by Drell. From that time, numerous attempts to model 
the form factor F 7r o 7 * 7 * and to get the predictions of the leading order decay rate within various 
approaches have been made [2, 3, 4, 5, 6, 7]. Recently this decay has attracted a renewed theo- 
retical interest in connection with a new precise branching ratio measurement by KTeV-E799-II 
experiment at Fermilab [8] with the result 

B(tt° -> e + e~(j),x D > 0.95) = (6.44 ± 0.25 ± 0.22) x 10~ 8 . (1.1) 



Here the Dalitz variable 



XD = Hd±- = l_2^L (1.2) 



(where E~ is the real photon energy) has been bounded from below in order to pick up the region 
where the final state radiation is soft and where the contribution of the Dalitz decay 7r° — > e + e~7 
which dominates at low xu is suppressed. Subsequent comparison with theoretical predictions of 
the SM based on the dispersive approach and various models for the pion transition form factor 
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(including the CELLO [9] and CLEO [10] data) has been done in [11]. The necessary ingredient 
of such an analysis is a good understanding of the QED radiative corrections to the process. The 
KTeV analysis used the early calculation of Bergstrom [12] to extrapolate the full radiative tail 
beyond xd > 0.95 and to scale the result by the overall radiative corrections to get the lowest 
order rate with the final state radiation removed with the result 

B™-*f(<!r° -> e+e") = (7.48 ± 0.29 ± 0.25) x HT 8 . (1.3) 

This should be compared with the SM theoretical prediction of [11] which has been found to be 
almost insensitive to the model dependent part within the relevant class of models for the form 
factor F„-o 7 * 7 *. Using the CLEO+OPE they obtained 

B^ iad (f° -> e + e~) = (6.23 ± 0.09) x 10~ 8 . (1.4) 

The result of the analysis can be interpreted as a 3.3a discrepancy between the theory and the 
experiment. This discrepancy initiated further theoretical investigation of its possible sources. 
Aside from the attempts to find the corresponding mechanism within the physics beyond the SM 
[13, 14, 15, 16] also the possible revision of the SM predictions has been taken into account. The 
theoretical estimate of the mass corrections to the decay width using the Mellin-Barnes represen- 
tation has been made in [17, 18] and this effect has been found to be negligible (the central value 
of the SM prediction is shifted by 0.5%). Also the incorporation of the new BABAR data [19] on 
the semi-off-shell form factor i ? 7r o 77 » in the time-like region into the analysis [20] has not influenced 
the SM prediction (1.4). 

The QED radiative corrections as a possible source of the discrepancy have been revisited 
calculating the contributions of the vertex-, box-type and self energy two-loop graphs in the dou- 
ble logarithm approximation [21]. The result has occasionally confirmed quantitatively the old 
Bergstrom calculation [12] which used a different type of approximation based on shrinking of the 
one-loop leading order graph into a local ir°e Jr e~ vertex. 

The aim of our paper is to present a more detailed analysis of the two-loop QED radia- 
tive corrections without using any approximation in order to check the validity of the previous 
approximative results. The natural formalism to treat the problem systematically is the Chiral 
perturbation theory (xPT) [22, 23, 24] enriched by photons and leptons [25, 26]. The leading order 
amplitude which is 0(u 2 p 2 ) within the chiral power counting has been calculated in [4] and the 
matching of the relevant low energy constant to the QCD in the leading order of the large Nq 
expansion has been done in [6]. The next-to- leading order contributions have not been calculated 
within this formalism yet. They can be divided into two groups. The first group corresponds to 
the additional strong higher order corrections to the 7r°7*7* vertex with pions inside the loops and 
it counts as 0{a 2 p A ) while the second one collects the pure QED corrections of the order 0(a 3 p 2 ). 
It is the latter group we will concentrate on in this paper. The relevant contributions consist 
of the six two-loop Feynman diagrams, namely the one box-type, two vertex-type, the electron 
self-energy insertion (these have been approximately investigated in [21]) and two vacuum polar- 
ization insertions. In order to renormalize the one-loop UV sub-divergences the corresponding 
one-loop counterterm diagrams have to be taken into account. Finally the remaining superficial 
UV divergence has to be renormalized by tree counterterm graph. The box-type diagram suffers 
further from the IR divergence, this is cancelled within the inclusive ir° — > e + e~{j) width. 
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In this paper we address the technical aspects of the calculation of the six two-loop Feynman 
diagram contributions mentioned above. The standard strategy consists of their reduction to the 
dimensionally regularized scalar integrals which will be subsequently expressed in terms of the 
eighteen Master Integrals. This can be done using the Laporta-Remiddi algorithm [27, 28] which 
is based on the integration by parts identities [29, 30] and Lorentz invariance identities [31]. We 
then calculate the Master Integrals using the technique of differential equations [32, 33, 34, 35, 36] 
and expand them up to and including the order 0{e) (where e = 2 — d/2) in terms of the harmonic 
polylogarithms [37]. Some of the Master Integrals has been already published in the existing 
literature, we either take them over [38, 39] or make independent calculations in alternative bases 
within individual topology classes and afterwards check the results [40, 41, 42, 43, 44, 45, 46]. This 
re-calculation found agreement with the formulae published earlier. We have also added new yet 
unpublished parts of some of the Master Integrals (typically the 0(e) terms of their e— expansion) 
in the closed form for the first time. We also discuss in detail the aspects of the UV renormalization 
of the two-loop graphs including the counterterm graphs described above and the treatment of 
the IR divergences within the soft-photon approximation. We give also the numerical analysis and 
discuss the various approximation to the exact two-loop expression. 

The paper is organized as follows. In Section 2 we summarize our notation and discuss the 
general structure of the amplitude. The third section is devoted to the general aspects of the 
systematic chiral expansion of the amplitude. Here we also give a list of the two-loop and one-loop 
Feynman diagrams contributing to the next-to-leading order pure QED corrections. In Section 4 
we discuss the general strategy of the renormalization of the one-loop and two-loop contributions 
and the treatment of IR divergences within dimensional regularization in detail. Section 5 is 
devoted to the calculation of the one-loop graphs with one one-loop counterterm vertex and also 
our renormalization scheme is specified there. In Section 6 we calculate the two-loop graphs and 
in Section 7 we discuss the soft-photon bremsstrahlung. In Section 8 we put all the ingredients 
together and give the final result for the virtual and real QED radiative corrections. We also 
discuss large logarithm approximation and relate our result to the Bergstrom's calculation. Some 
preliminary phenomenological applications are discussed in Section 9. In Section 10 we give a 
brief summary and conclusion. Some technical details are postponed to the Appendices. The 
relevant part of the xPT Lagrangian with virtual photons and leptons is summarized in Appendix 
A. The reduction of the six two-loop graphs to the scalar integrals is presented in Appendix B. 
In Appendix C we list the integration-by-parts identities for the scalar integrals and in Appendix 
D we summarize the results of our (re-)calculation of the relevant Master Integrals and give a 
comparison with existing literature. 

2 Basic properties of the amplitude 

In this section we discuss the basic features of the amplitude. We set the notation and kinematics 
and then briefly comment on the general properties of the lowest order amplitude which corresponds 
to 0(a 2 ) order in electromagnetic interaction (and all orders in QCD for the pion transition form 
factor, which is here the only nonperturbative ingredient). 
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2.1 Notation and kinematics 

The invariant amplitude -M^o^e+e- for the decay is defined by means of the matrix element 

(e+(q + , «+)e-(g_, out|7r°(Q);in> = i(2vr)V 4) ((5 - q+ - qJ)M^ e+e - (2.1) 

which is supposed to be calculated in the presence of strong and electromagnetic interactions. 
According to the Lorentz covariance we can further write 

M n o^ e+e - =u(q-, s-)T n o e+e -(q-,q + )v(q-, s_) (2.2) 

where r n o e + e -(q-,q+) is a one particle irreducible 7r°e + e _ vertex. Off shell it can be conveniently 
decomposed introducing four scalar form factors P, A± and T defined as 1 



ii> e+e - (q.,q + ) = P(q 2 _,q 2 +, QV + (0_ " m)^ A.(q 2 _,q 2 + , Q 2 ) 

+ A + (q 2 _,q 2 + , Q 2 )7 5 (^ + + m) + T(q 2 _,q 2 + , Q 2 )(0_ ~ m) 7 5 (^ + + m). 

Here Q = q+ + q- and the charge conjugation invariance implies 



(2.3) 



A_(q 2 _,q 2 + ,Q 2 ) = -A + (q 2 + ,q 2 _,Q 2 ) 
P(q 2 _,ql,Q 2 ) = P(q 2 + ,qtQ 2 ) 

T(ql,q 2 + ,Q 2 ) =T(q 2 + ,q 2 _,Q 2 ). (2.4) 

For the electron-positron pair on shell we get then 

iM n o^ e+e - =u(q-,S-)~f 5 v(q^, S -)P(m 2 ,m 2 ,Q 2 ), (2.5) 

and, as a consequence, the total decay rate is given solely in terms of the on-shell form factor 
P(m 2 ,m 2 ,M 2 ) as 



r>_> e+e - = ^P(M%>) \P(m 2 ,m 2 ,M 2 )\ 2 (2.6) 



where 



is the velocity of the electron-positron pair in the CM frame. 

Note that the semi-on-shell form factor P(m 2 , m 2 , Q 2 ) can be extracted from the one particle 
irreducible vertex T n o e + e -(q-,q + ) by means of the following projection 

P(m^m^Q 2 ) = - lim ^Tr[(^_ + m)r 7r o e+e -(g_,g + )^ + -m)7 5 l . (2.8) 



i±- 



This dimensionless form factor is an analytical function of the variable s = Q 2 in the complex 
plain with a cut [0, oo) where the unphysical threshold Q 2 = corresponds to the two-photon 
intermediate state. For further convenience we introduce two dimensionless kinematical variables, 
namely 

Q 2 

y 4m 2 ^ ^ 



1 In what follows we use the convention e 0123 = 1 and 75 = i7°7 1 7 2 7 3 
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Figure 1: The leading order 0(a 2 ) contribution to the amplitude. The blob corresponds to the 
pion transition form factor. 



and 



l-i-i 



I = W)_i = v ^ (2 , 0) 

which map the unphysical threshold to y = and x = 1 respectively. 

In what follows we assume perturbative expansion of the amplitude in the QED coupling a. 
Consequently we can write for the form factor P 

P(m 2 ,m 2 , M 2 ) = P LO (m 2 ,m 2 , M 2 ) + P NLO (m 2 , m 2 , M 2 ) + 0(a 4 ), (2.11) 

where P LO = 0(a 2 ) and P NLO = 0(a 3 ), and for the decay rate 

r(7T° -> e+e") = r LO (7r° -»• e+e") + T nlo (tt° -»• e + e~) + O(o 6 ). (2.12) 

In order to cancel the infra red (IR) divergences present in Y NLO we have to add also the real photon 
bremsstrahlung contribution and consider inclusive decay rate of the process 2 7r° — > e + e~( r y). The 
size of the NLO real and virtual QED radiative corrections can be then described by means of 
the factor 5(x c ] ^ t ) defined as 

r NLO {7r° -> e + e-{ 7 ),x D > x^) = 5(a;^ ut )r iO (7r -> e+e") (2.13) 

where x/j is the Dalitz variable (1.2) and where all the tt° — > e + e _ (7) events with xn > xff are 
included. 



2.2 The amplitude at the order 0(a 2 ) 

Within the QED the leading order contribution P LO (m 2 ,m 2 ,M 2 ) to P(m 2 ,m 2 , M 2 ) is of the 
order 0(a 2 ) and comes from the diagram depicted in Fig. 1. The bubble there corresponds to the 
strong matrix element 

-e 2 J d 4 x e ik - x (0\T(f(x)f(0))\n°(Q)) = ie 2 s^k a Q p F 7r o rr (k 2 , (Q - k) 2 ) (2.14) 

2 Note that the same final state has also the Dalitz decay tv° — > 77* — > ye + e~, which is however dominant in 
different region of the phase space corresponding to small xd (cf. (1.2)). For large enough id the Dalitz decay 
contribution is tiny, however, for the experimentally used cut oniD this contribution should be also included. 
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where j^{x) is the hadronic part of the electromagnetic current 



f = h Sr f u -'^drfd (2.15) 



and i ? 7r o 7 * 7 * (k 2 , I 2 ) is the pion transition form factor. 
The explicit formula reads then 



G + g+)aG-g-)/3 1 rglfi v 

X ((I - g_) 2 + i0 )[(Z + g + ) 2 + iO] 7 ^ / - m + iO lv ' 1 &j 
and using the projection (2.8) we get (cf. [47]) 

P LO (m 2 m 2 O 2 ) - -i— [ dH F ^*r( D{ ~ ] > Di+) ) x(0 2 D (-) D ( +)) , 2 m 



where we abbreviated 



D (±) = (I ± g±) 2 + iO, (2.18) 
£>(°) = Z 2 - m 2 + iO (2.19) 

and where 

A(a, 6, c) =a 2 + b 2 + c 2 - 2ab - 2ac - 26c (2.20) 

is the triangle function. 

The pion transition form factor represents the unknown nonperturbative QCD ingredient of 
the above formula, therefore the evaluation of the integral (2.17) is model dependent. However 
some general features of (2.17) can be deduced in a model independent way. We will discuss them 
in the rest of this section. 

Let us first briefly remind the properties of i ? 7r o 7 * 7 *. It has the following short distance 
asymptotics which is a consequence of OPE. For Q fixed and A — > oo we have [48] (see also [7]) 

F n o rr ((Xl) 2 , (Q - XI) 2 ) = - prW^r (1 + 0(a s , A" 1 )) (2.21) 

(where F n is the pion decay constant) and therefore we can conclude that the integrals (2.16) and 
(2.17) are convergent. Also the long distance asymptotics of F n o~*~*(l , (Q — I) 2 ) is known from 
the first principles being fixed by the QCD chiral anomaly. Namely in the chiral limit 

^(0,0) = ^-, (2.22) 

and therefore the low energy behavior is expected to be given by the chiral expansion of the form 

F^. ri P, (Q - if) = ^ (l + O (£, . (2.23) 

where Ah ~ lGeV is the hadronic scale limiting the applicability of \PT. 

Another useful model independent property is that the two-photon intermediate state contri- 
bution to the imaginary part of the leading order form factor (2.17) with Q 2 = s extended off the 
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pion mass shell 3 is uniquely fixed by the value F 7r o 7 * 7 »(0, 0). Indeed, using the Cutkosky rules and 
cutting the two internal photon lines we get [1],[2] 



Im P LO (m 2 ,m 2 ,s)\ 2l = 9(s)ira z m- 



l-fi(s) 

0{8) ~\1 + P(8) 



In 



i^o 7 * 7 *(0,0). 



(2.24) 



For s < (M n o + 2M n +) 2 this contribution saturates the imaginary part of P (m 2 ,m 2 ,s). As a 
consequence, the known imaginary part Im P LO (m 2 , m 2 , M 2 ) and the known width F^P of the 
27 pion decay in the leading order in the QED expansion, 



t 7 = ivrQ 2 M3 () |F w o rr (0,0)| 2 , 



can be further used to get another exact result concerning the branching ratio 



B lo (tt° 



R 



B LO (7T° -»• 77) 



(2.25) 



(2.26) 



Namely, using InLP i0 | 27 instead of P LO in (2.6), we get the following model independent unitarity 
bound [21 



R > 



1 / am 



1 



2 \M w0 J /3(M 2 )"' \l + P(M 2 ) 



In 2 



4.75 x 10" 8 . 



(2.27) 



The explicit knowledge of Im -P LO | 27 allows also to pinpoint the most important nonanalytic 
contribution to P LO (m 2 ,m 2 ,s), namely that stemming from the two-photon intermediate state. 
The latter is given by means of the following once subtracted dispersive representation [3] 



t^to, 2 2m n m, 2 2 «m s f°° ds' Im P LO (m 2 ,m 2 ,s' ) 1 27 
P i0 (m 2 ,m 2 ,s)| 27 = P LO (m 2 ,m 2 ,0)| 27 + - / — !_>JJI, 

TJO S S - S 

and it is therefore fixed uniquely up to one unknown subtraction constant P (m 2 , m 2 , 
dispersion integral for the physically relevant region s = Q 2 > Am 2 reads 



P LO (m 2 } m\ s )\ 2 disp = 9l r 

7T Jo 



ds Im P LO (m 2 , m 2 , s' ) 1 2 7 



a 2 mF 7r o 7 * 7 *(0, 0) 



s - Q 2 - i0 
1 



•0(8) 



Li 2 (x) - Li 2 ( - ) + m In (-x) 



(2.28) 
| 27 . The 

(2.29) 



Note that _F 7r o 7 * 7 » (k ,1 ) can be obtained by means of the LSZ reduction formula from the three-point correlator 



where 



d 4 x d i ye ik ^ Q - k ^(0\T(j^x)f(y)P 3 (0))\0} = k a Q?H(k 2 , (Q - k) 2 ,Q 2 ) 



Pa — \ a 5 

= 7 9 



is the pseudoscalar density. Namely for Q 2 — > M 2 we have 

H{k?,{Q-k)>,Cf) = |iZ^^jM^ 07 , 7 ,( fc 2 ,(Q- fc f) 
+0 ((Q 2 - M 2 )°) . 

This offers the natural possibility to extend the form factor F n o J * J * (k 2 , 1 2 ) off the mass shell 



l2 T\yf2 



p (h 2 (n — b\ 2 n 2 \ — 2 (Q ~ M^q) 

^o 7 . 7 .(« ,w K i >w ; - 3 ( ff Q(Q)|p3( )|o) 



H(fc 2 ,(Q-fc) 2 ,Q 2 ) 
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In this formula x < is given by (2.10) and LI2 is the dilogarithm defined as 



Li 2 (z) 



I jln(l-t) (2.30) 



which is analytic in the complex plain with cut [l,oo). The unknown explicit form of the form 
factor i ? 7r o 7 . 7 * in the intermediate energy region can influence only the subtraction constant 
P LO (m 2 ,m 2 ,0)| 27 . 

As a consequence, we can split the leading order amplitude P LO (m 2 , m 2 , s) into two parts, 
namely 



P LO (m 2 ,m 2 ,s) = P LO (m 2 ,m 2 ,s)\ 2l4isp 
+2a 2 mF 7r o 7 * 7 *(0,0) 



3 . / ?7i 2 \ 5 { s m 2 



2 ln [ A 2 J 2 + X \ A 2 ' A 2 



(2.31) 



The first part corresponding to the two photon intermediate state is completely independent on 
the details of the pion transition formfactor. The second part (where A is an intrinsic scale 
characteristic for the formfactor i ? 7r o 7 » 7 », i.e. the scale at which is the integral (2.17) effectively 
cut off) accumulates the above subtraction constant P LO (m 2 ,m 2 , 0)|2 7 as well as the contributions 
of higher intermediate states which appear when also the bubble in the Fig. 1 is cut. The explicit 
form of the x{~fi7i jh) as a functional of unknown transition formfactor i ? 7r o 7 * 7 * has been derived 
in [18] with use of the Mellin-Barnes representation. 

We have explicitly kept apart the logarithmic term in the square brackets of (2.31). The reason 
is that this term corresponds to the leading dependence of P LO (m 2 ,m 2 , s) on the effective cut-off 
scale A. The origin of this term is easy to understand [3]. In the case of point-like pion (i.e. when 
F 7r o 7 » 7 » (k 2 , 1 2 ) = F 7r o 7 » 7 » (0, 0) = const.), the integral (2.17) is logarithmically divergent. Using 
the sharp cut-off at the scale A instead of effective cut-off provided by i ? 7r o 7 * 7 * (I 2 , (Q — I) 2 ) we get 
for A — > 00 

P LO (m 2 ,m 2 ,s)\fZV^° S = -3a 2 mi^ r (0,0)ln (^) + O(l) (2-32) 

where the O(l) part includes terms that are finite or suppressed for A — > 00. Thus we expect the 
same behavior also for the full amplitude P LO and therefore x = 0(1) for A — > 00. 

Because the discontinuities of x as a function of s start at s ~ A 2 , \ is analytic in the physical 
region s < A 2 and can be therefore expanded in this region in the power series of the variable 
s/A 2 . This suggests that x(j7,^r) 

can be approximated at the leading order of this expansion 

by A- independent constant. 



3 Systematic chiral expansion 

Within the SU(2) x SU(2) variant of xPT supplemented with dynamical photons and electrons 
(cf . [26] ) the formfactor P is given in terms of the systematic simultaneous expansion in powers of 
the momenta (s), the quark and electron masses and the fine structure constant a, to which the 
chiral orders are formally assigned according to 

s, m q ,m 2 ,a = 0(p 2 ). (3.1) 
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Figure 2: The leading order 0(a 2 p 2 ) contribution to the amplitude. The shaded box corresponds 
to the ne + e~ counterterm contribution. 

The relevant parts of the enlarged X-PT Lagrangian are listed in the Appendix A. The hierarchy 
of the various contributions is then controlled by the Weinberg power-counting formula [22]. As 
is usual in xPT, for the regularization of UV as well as IR divergences we use the dimensional 
regularization (DR) in what follows. In order to avoid problems with intrinsically four-dimensional 
objects like Levi-Civita pseudo-tensor and 75 we use here the variant known as Dimensional Reduc- 
tion. For calculation of P this means that we first project out the form factor from the amplitude 
by means of (2.8) using four-dimensional Dirac algebra and only then we dimensionally regularize 
the resulting scalar integrals. 

3.1 The leading order of the chiral expansion 

The leading order contribution P xLO (m 2 , m 2 , s) represents the chiral order 4 0(a 2 p) and can be 
divided into two parts (see Fig. 2). The first one is a logarithmically divergent one-loop graph 
with local 7r°7*7* vertex stemming form the 0(ap 2 ) Wess-Zumino-Witten Lagrangian (here and 
in what follows we write down only the relevant vertices) 

*tfzw = I O Y^ aP F^F^ + ... (3.2) 

(Fo is the pion decay constant in the chiral limit) and the second one corresponds to a tree-level 
counterterm graph originating in the 0(a 2 p 2 ) Lagrangian 



'-'■nee ~ H 1 



7T 



X r (/-0+o (-+ln4vr 



2 V £ 



<9„7T° 

e^^e^— + • • • , (3.3) 



where x r (^) is a renormalized counterterm coupling at a scale fi. One can think of the loop part 
of P xLO as approximating the leading order formfactor P LO given by (2.17) by means of inserting 
the leading order term of the chiral expansion of the pion transition formfactor (2.23) into (2.17). 
This insertion however modifies significantly the high energy region of the loop integration starting 
at the onset of the resonances where the chiral expansion fails to converge. Such a modification of 
the loop has to be compensated by local counterterm contribution in such a way that the 0(a 2 p) 
term of the chiral expansion of the P LO is exactly reproduced. This is the general idea of the 



4 Note, that for the amplitude M jr o^ e + e - we have M n o^ e + e - = 0(a 2 p 2 ), because the fermion wave functions 
are counted as 0(p 1 ^ 2 ). 
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matching of the coupling constant x'X/- 4 )- I n the absence of the first principle determination of 
F 7r o 7 * 7 » this matching procedure is however model dependent. We use here the value 

X r {^ = 770MeV) = 2.2 ± 0.9 (3.4) 

which has been obtained in [6] by means of using a large Nc inspired Lowest Meson Dominance 
(LMD) ansatz for F n o~*~*. 
As a result we get 5 



/ a \ 2 2m 



3, fm 2 \ 5 



(3.5) 



The structure of this expression can be easily understood. It corresponds to the formula (2.31) 
where both F 7r o 7 . 7 » (0, 0) and x(s/A 2 , ?n 2 /A 2 ) have been replaced with the leading terms of their 
chiral expansion, provided we identify the renormalization scale /i with the intrinsic cut-off scale 
A of the formfactor F^o^*^* . 

3.2 The 0(a 3 p 2 ) part of the next-to-leading order 

At the next-to- leading order the amplitude is generically 0(p 8 ) in terms of the simultaneous 
expansion according to the chiral order assignment (3.1). At this order there are two types of 
contributions, which counts either as 0(a 2 p 4 ) or as 0(a 3 p 2 ) for the amplitude (or 0(a 2 p 3 ) and 
0(a 3 p) for the form factor P). The contributions of the first type collect two- loop and one- loop 
graphs with the same topology as depicted in the Fig. 1 (where now the blob represents either 
a one loop subgraph with pion internal lines or 0(ap 4 ) order counterterm) and in addition tree 
graphs with 0(a 2 p 4 ) counterterms 6 . Such contributions can be understood as was mentioned 
above as the next-to-leading terms of the chiral expansion of (2.31). 

In this paper we will concentrate on the contributions of the second type (i.e 0(a 3 p 2 )) which 
represent the pure QED corrections. In this case we get again three classes of graphs, namely 
six two- loops graphs with virtual photons and electrons (see Fig. 3), six one-loop graphs with 
0{ap 2 ) counterterms (see Fig. 4, graphs (2)-(6)) or 0(a 2 p 2 ) counterterm (see Fig. 4, graph (1)) 
which renormalize the one-loop subdivergences of the corresponding two-loop graphs and 0(a 3 p 2 ) 
tree- level graphs which are necessary to renormalize the remaining superficial divergences. The 
latter are of the same order in p as the counterterm (3.3) and therefore the relevant vertex from 
the 0(a 3 p 2 ) La grangian can be summarily written in the form 

Cf£ = (f) 3 + 0(e~ 2 ) + O^- 1 )] ef^e^, (3.6) 

where £ r (/x) is the renormalized coupling and we have not written the UV divergent part explicitly. 



5 At this order we can put Fa = F v . In fact, such a replacement corresponds to partial re-summation of the higher 
order corrections, namely the renormalization of the pion external leg. 

6 Strictly speaking there is also contribution from the pion external leg renormalization, which we have however 
effectively added in the LO by means of the replacement Fq — > F^. 



11 




Figure 3: The two-loop part of the next-to- leading order 0(a 3 p 2 ) contributions to the amplitude. 

4 Structure of the two-loop corrections and renormalization 

In this section we briefly discuss the general structure of the one and two-loop contributions within 
the dimensional regularization and within renormalization scheme suitable for power counting non- 
renormalizable effective field theories [49] (cf. also [50]). This discussion will be helpful for the 
organization of the results of the explicit calculation and for consistency checks of the results 
presented in the next sections. 

Let us write the effective Lagrangian in the form 



where C n accumulates the counterterms which are needed in order to renormalize the superficial 
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(1) 



(2) 



(3) 




Figure 4: The one-loop part of the next-to-leading order 0(a 3 p 2 ) contributions to the amplitude. 
The numbers labeling the graphs are in one-to-one correspondence with Fig. 3. The shaded box 
is the 0(a 2 p 2 ) counterterm while the blobs represent 0(ap 2 ) ones. 

divergences of the n— loop graphs and fx is the dimensional regularization scale. Schematically 

0«, (4.2) 

where i^^(//) n_loop are the renormalized counterterm couplings (i.e. the finite parts of the coun- 
terterms at the n— loop level) and 0$ is a set of operators. In the above formula for C the 
factor n~ 2en naturally appears when the scale [i is artificially introduced into each loop integration 
writing d d k = fi~ 2£ (^ 2£ d d k) in order to restore the four-dimensional canonical dimension of the 
d— dimensional integration measure. 



#r^n-loop 



n 

£ 



.,(0. 



-loop 



4.1 Renormalization of the one-loop contributions 



For the total one-loop contributions to any one-particle irreducible vertex 7 entering the game we 
can write schematically 



1 



1— loop 



yT 2£ m D 



2e 



' 1— loop 

7-i 



+ 7 



1— loop 



+ £ 7i" l0 ° P + 0(^) 



/ l-loop ' 

+ \^)-'h±— 



(4.3) 



where D is the (four-dimensional) canonical dimension of 7. In this formula the first line represents 
the result of the loops while the second line accumulates the counterterm contributions. In this no- 
tation, both 7i7 loop and x^i^) are generally polynomials in the external momenta. Moreover X^G") 
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is linear in the renormalized counterterm couplings K^(/j) 1 loop introduced above. On the other 
hand, the functions 7 1 loop for i > — 1 have more complicated analytical structure with branch 



0— loop 



points and cuts. In general, j- ~ oop depend nonlinearly on the tree-level couplings KT(fi) 
The chosen normalization ensures that 7 J 1_loop are dimensionless. The above structure (4.3) is 
shared by the LO contribution to the formfactor P and also by the UV divergent one-loop sub- 
graphs of the NLO corrections to P, namely the off-shell one- loop ire + e~ and je + e~ vertices, the 
electron self-energy and the vacuum polarization. 

We use here the renormalization scheme suitable for power-counting nonrenormalizable ef- 
fective theories [49] and require renormalization scale independence order by order in the loop 
expansion. For -y 1- loo P this means that the following finite quantity 



Xt = — 



-2c- 



1— loop 1 

7-i 



+ 



1— loop 

7-i 



X»+7^° OP ln(^) 



(4.4) 



has to be ^/—independent. This implies a running of x^in) according to 



- 7 ii l0 ° P In ( £2 



+ e 



X 7 In 



m- 



1 l-loop i 2 

o7-i In 



m. 



+ 0(e 2 



(4.5) 



and corresponding running of the (linear combinations of ) counterterm couplings K^(fi) 1 loop . As 
a result, (4.3) can be re-organized in the following simple manifestly RG invariant form 



7 l-loop =m D 0-loop + _^ + 0(£; 



(4.6) 



4.2 Renormalization of the two-loop contributions 

For simplicity let us first assume that there are no infrared (IR) divergences. The generalization 
with the presence of IR divergences will be discussed in the next subsection. 
At the two loop level, we have three types of contributions, namely 



2-loop 2-loop . l-loop . tree 

7 = 7l + Ict + 7cr • 



(4.7) 



The first one corresponds the genuine two-loop contribution which can be schematically written 
in the form 

Ar _ /„,2— loop ^ 2— loop 

2— loop 



2— loop 

7l 



. / 2— loop 2— loop 



(4.8) 



The second one represents a sum of one-loop graphs with one-loop level counterterms which are 
necessary to renormalize the sub divergences of the two-loop part 



l-loop -4e D ( M 

ICT = V \ 



2c 



(i), l-loop- 
l-loop _ '-1 



/ ^(i), l-loop 

x + C®> 1-10015 + eC?' 1 " loop + 0(e 2 ) 



(4.9) 



where we have explicitly pulled out the dependence on the renormalized one-loop counterterm 



couplings x^(fi) 1 loop and coefficient of the corresponding infinite parts 7_{ 



(i), 1— loop 



Finally we 
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have also a tree level contribution of the counterterms necessary to renormalize the remaining 
superficial divergence of -y 2_lo °P + ^ T loop 5 naively 



tree _ -4e D 
1CT — f 1 m 



CM 



2— loop ^ (*)> 1— l°°P/^((*): 1— loop 
7-2 _ Z^j7-i 



£- 



(i), 1— loop^-,(i), 1— loop 







<(^) 1 " 1 ° OP G-{ 



(i), 1— loop 



9 2-loop ^ (i), i-ioop^i;, 

z 7-2 Z^i 7-i 



(i), 1— loop^-y(i), 1— loop 



In 



(4.10) 



(here ^(//) is polynomial in external momenta and linear combination of the couplings x[(/i) 2_loop 
from £2)- However, in the absence of the IR divergences, in the sum 7^ loop + 7( ^ r loop onry 
the local (i.e. polynomial in the masses and external momenta) UV divergences survive. This 
fact implies nontrivial relations between various 7? -loop ) cj 1 loop and 7^' 1 loop which can be 
used either in order to simplify the above contributions or as a nontrivial check of the explicit 
calculations. Cancellation of the explicitly \x— independent nonlocal 0(e~ 1 ) terms in the sum 

2— loop . 1— loop 

7l 

divergence of the one-loop graph) 

7_i 



+ 7ct°° p needs (note that C^\' 1 loop has to be local because it corresponds to the UV 



2— loop V""^ s~t(i), 1— loop (i), 1— loop 



7ll 



2— loop 

7-i 



(4.11) 



where y^i° op j is local. Similarly, cancellation of the nonlocal 0(s 1 ) terms proportional to 
In (/i 2 /m 2 ) requires 

27 2-ioo P _ 7 W> i-k^W- ^ 1oo p = o. (4. 12 ) 



In fact these relations are valid also graph by graph. Introducing RG invariant counterterm 
couplings Xi according to (4.4), namely 



1— loop 



-2e 



(i), 1-loop' 
1-loop _ 2^1 



+ 



(i), 1— loop 

7-i 



(4.13) 



we can write using (4.11) and (4.12) 



2-loop . 1— loop — 4 e D ( t 1 

7i + Ict = V m — 



lc 



2-ioo P c { i ] { 1 - lo °p- i - l °°p 



7-2 



2— loop 

+ 7-i 



+ 



+ 0(e) 



(4.14) 



and the remaining tree-level contribution coming from the two-loop counterterm can be then 
simplified as 

2- 



7 ^ e = ^"V (ii 



l£ 



2-loop £ C { 1 ] { 1 - loop x j 1 - lo ° P 



7-2 



ct ov ) i 



(4.15) 
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Here the UV finite and RG invariant ' local quantity f „ is a two-loop analog of x 



-4e 



2 loop ( 7 -l°°P) + £. x r 1-loop 



7 2-ioo P ( 7 ^ ioop ) _ + y\ ^- ioop cW' 1-10015 



+ 



if! 

m 2 



2 ( 7 !- loop ) + ^<( /U ) 1 - looD ci i J' 1 ~ loop 



+7 - loop ln^£ f j + 0(,). 
This implies the following explicit dependence of on the renormalization scale \x: 



2-loo Pl 2 ( A* 



i 

2 



1 _1— loop^-,(i), 1— loop 



In 



(4.16) 



(4.17) 



As a result we get the following simple and manifestly RG invariant form of the total two-loop 
contribution 7 2 ~ 1o °p to 7: 



^2-loop _ m D 



To 



2— loop 



-,(«), 1— loop_l— loop /^(*), 1— loop (i), 1— loop\ 



(4.18) 



4.3 Treatment of IR divergences 

The presence of the IR divergences complicates the above simple picture a bit. After the UV di- 
vergences are subtracted, additional divergent terms survive, generally both in 7^ _loop and 7ct°° p , 
namely 



/ 2-loop 2-loop ' 



(4.19) 



and 



(i), 1— loop 



7^ = M-v(|fE c 



£ 



(i), 1-loop' 
1-loop _ J^l 



, /nr(i)i 1-loop / (j) 1-loop' 



(4.20) 



However, in the sum of these two contributions, the 0(e 2 ) part has to vanish. This gives another 
useful relation which can be used as nontrivial check of the explicit calculation, namely 



2— loop V~~* x-i(i): 1— loop (i), 1— loop 

~ / ^-l.IR 7-1 



(4.21) 



7 As in the previous subsection we require RG scale invariance order by order in the loop expansion. We also tacitly 
assume that all the relevant contributions are included in 7 which represents a RG invariant physical observable. 
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and we get for the remaining IR divergent part 



1lR 



fx m — 



2 — loop . V"^ /"<(!■) 

2-loop _ -4s D ( M \ 6 + °-l 



1— loop_l— loop 



,IR 



(4.22) 



In our case the only one one-loop counterterm graph with IR divergence is that corresponding 
to the counterterm of the 7re + e~ vertex graph (let us denote the corresponding C^^ R loop as 



C^l' IR loop in what follows) and in (4.22) the only relevant x\ loop we denote as x- The IR 
divergences in \^/\ 2 where 



— 1— loop 



171 = 17 



l-loop|2 _|_ 2-loop f 1-loop \ _|_ 1-loop / 2-loop 



(4.23) 



can be cancelled to the given order including the soft bremsstrahlung contribution Abs into the 
inclusive decay rate, where 

Abs = |7 1_loop | 2 ^s (4.24) 

and where Ibs is a dimensionally regularized phase space bremsstrahlung integral (see Section 7) 
with a general structure 

Ibs = ^ + h + 0(e)) . (4.25) 

It holds schematically 



7/« l00P (7^ loop )* + 7^ loop (jm°° P Y + |7 1_loop | 2 ^s = 0(e c 



(4.26) 



Inserting (4.6) and (4.22) into (4.26) and collecting the coefficients at various powers of \ in the 
0{e~ 1 ) term gives the following conditions 

J_i + 2Re(ci 1 1 ) / ^ loop ) = (4.27) 



and 



Re^ 1 - 1 ^)/-! +Re(7o^ loop C^; / ^ luop ^) + Re( 7 l7; oop ) = 



1— loop^l), 1— loop*N 



2-loop n 



(4.28) 



|7o~ lo0 T/-i + 2Re( 7o 1 - loop 7_T loo f ) = 0. 



The latter two relations can be rewritten with help of (4.27) as 

2— loop /o(l)> 1— loop 1— loop 

7-1, IR — U -1,IR 7o 



(4.29) 



(4.30) 



which represents another relation which can be used as a check of the explicit results. Using these 
relations, we get finally 

2 



1— loop i — 

7o +X 



(1 + J ) - 2Re 



1— loop . — \ 1— loop ,^-,(1), 1— loop 

7o +X 7i l C y ' R 



+2Re<! (7o~ 100P + X 



2-loop -p ( n {i), 1-loop 1-loop n {i), 1-loop (£), 1-loop 

7o + s + 2-^i V o x i ~ °i 7-i 



(4.31) 



which is manifestly independent on the RG scale \x. 
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5 The one-loop graphs 



In this section we summarize the results for the one-loop graphs which are relevant for the full 
0(a 3 p 2 ) two- loop calculation. First we will present the one-loop contribution to the form factor 
P(s, m 2 ,m 2 ) in more detail including also the order 0(e) which will be necessary for treating the IR 
divergences. We will also discuss the one-loop UV divergent subgraphs of the two-loop graphs and 
one-loop graphs with one counterterm vertex. We will explicitly point out the individual orders 
in e according to the general structure discussed in the previous section. In order to simplify the 
results and to avoid some repeating multiplicative factors, we will present all the results in terms 
of the re-scaled form factor 7(2) defined as 



j(s,m 2 ,m 2 ) = 2— (-Y P(s,m 2 ,m 2 ). (5.1) 



5.1 The leading order amplitude revisited 

As we have seen from the general formula (4.31), we need more detailed information on e— ex- 
pansion of the LO amplitude ■y 1 ~ loo P (namely the 0(e) term). In order to be consistent with the 
two-loop calculations it is also convenient to rewrite ^y 1 ^ 00 ? in terms of the Harmonic Polyloga- 
rithms of Remiddi and Vermaseren [37]. Let us note that in the physical region the kinematical 
variable x is negative (cf. (2.10)), while the two-loop integrals are originally calculated in their 
analyticity region corresponding to < x < 1 and only then analytically continued (cf. Appendix 
D ). This continuation apart from generating imaginary parts brings about also additional minus 
signs into the arguments of Harmonic Polylogarithms and therefore the most convenient way how 
to present the result in the physical region is to use as an argument of these functions a new 
variable 

z= - I= Trf >a (5 - 2) 

In order to further simplify the long expressions we use the notation 

7 = 7 — ln4-7r (5.3) 

and rewrite the rational function multiplying the Harmonic Polylogarithms H(a±, . . . ,a n ; z) in 
terms of the variable f3. 

As a result we get then for the one-loop contribution 7 1 ~ loop (z) (cf. Section 4 ) 



7 i-i°°P( 2 ) = ^ (£) ( 7o 1 - loop (z) + x + eil- l °° p (z) + 0(e 2 )) . (5.4) 
In the above formula (5.4) we have for < z < 1 

7d~ loop = 0; z ^ + i7rH ^ z ^ - 2H ^ ^ + t) - 1 + b (5 - 5) 

and the 0(e) term is 



7i 



1— loop 



7T 



7T 



2^ 



- ( y#(l; z) + -jH{-l; z) + 2F(-3; z) + —#(0; z 
-4#(1, -2; z) - 4ff(-l, -2; z) - 4ff(-2, -1; z) 
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-2H(-2, 0; z) + 2H(1, 0, 0; z) + 2H(-1, 0, 0; z) - H(0, 0, 0; z) 



ivr 2#(-2; z) + 2#(1, 0; z) + 2H(-1, 0; z) - H(0, 0; z) + 



7T 



7 2ff(-2; z) - iJ(0, 0; z) - mH(0; z) 



7T 



5C(3) 



7T 



+57 - -7 - — - 9 



where 



X 



'WIG 



2e X 2 7 2 U 2 



+ 0(e) 



For further convenience let us also mention explicitly the UV divergent part 



1— loop 

7-i 



The formula (3.5) can be then easily reconstructed as 



[m , m , s) 



2 \ttJ Rr e^0 ' 



(5.6) 



(5.7) 



(5.8) 



(5.9) 



and we can interpret the renormalization scale independent constant x m terms of the renormalized 
constant x r (/■*) a t scale \x = m as 



X = X r (m) - -7, 



numerically 



X 



■16.8 ±0.9. 



(5.10) 
(5.11) 



5.2 The one loop counterterms 

In this subsection we summarize the counterterms needed for renormalization of the one-loop sub- 
divergences of the two- loop graphs. We can write the relevant counterterm Lagrangian in the 
general form either in terms of the renormalized couplings (finite parts of the counterterms) 



£ = fi 



-2e 



71"/ 



- 1) verfD^e + (V 7 (/x) + ^ mee + Ufa) + ^\ ^i^, 



(here D = d + \eA) or in terms of the renormalization scale invariant constants (cf. (4.4)) 



c = „-* m m 



2e 



ttJ \mJ 



x 6 --) verfD^e +(x 7 + i) mee + fx 8 + — J -F^F^ 



(5.12) 



(5.13) 



These counterterms contribute to the electron self-energy S(p), the vertex function T{p,p) and 
vacuum polarization I4(p) which are related to the physical electron mass and physical charge, 
namely 



m 



phys 771 + Xl(771p^y S ) 



Cphys 



1 + 2 n (°) 



(5.14) 
(5.15) 
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Within our regularization scheme and at one-loop level we get 



S(m) 

n(o) 



\7T/ \m 



2c 



"4 7 "3 



7T 



771/ 



1 



:7 + ^8 



therefore adjusting 



4 1^-3 



x 6 - x 7 



(5.16) 
(5.17) 

(5.18) 
(5.19) 



(i.e. choosing the on mass shell renormalization scheme) we ensure, that the original parameters 
m and e coincide with the physical mass and charge. The parameter xq is connected with the 
electron wave function renormalization, namely 

0E(p) 



1 



(5.20) 



where in our scheme at one loop (here the e pole corresponds to the IR divergence) 



3E(p) 



a 

-J A* 

7T 



-2c / A* 



2c 



1 3 
'2^ + 4 



7 "3 



x 6 



(5.21) 



Therefore the parameter xq is not physical and has to cancel in the physical amplitudes. We can 
conveniently set 

51 / K\ 

(5.22) 



x 6 



7 



4 V 3, 

then the only effect of the electron wave function renormalization is the IR pole which is cancelled 
by the analogous pole in the diagonal part of the bremsstrahlung integrals (see Section 7). We can 
therefore forget the electron wave function renormalization completely provided we simultaneously 
throw away the IR divergent part of the diagonal bremsstrahlung integrals. 



5.3 One-loop graphs with counterterms 

There are six types of the one-loop graphs with 0(ap 2 ) and 0(a 2 p 2 ) counterterms which are 
depicted as (1) — (6) in Fig. 4. Let us denote their individual contributions as 7^ 1 loop . As a 
consequence of the symmetries of the graphs we have the relations 



(2), 1-loop _ (4), 1-loop (D, 1-loop 

<CT ~~ 'CT — 'CT 

(5), 1-loop _ (6), 1-loop (II), 1-loop 

/CT ~~ 7CT — 7cT 



(5.23) 
(5.24) 



and as a result of straightforward algebra we find that 7c T loop ^ and "yQ^°° P ^ are simply related 
to 7 1 ~ 1oo p ( c f. (5.4) and (5.8)), namely 




1— loop 

7-i 
e 



+ 7 



1-loop _,_ ^1-loop _|_ q^ £ 2\ 



+ £7i 



(n), i-ioop 



7ct 



(I 



-4e ( ^ 



■le- 



ft 



mJ \7T/ \ OS 




1— loop 

7—1 1— loop 1 1— loop 

— - — + 7 + £7i 



+ 0(e 2 
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(5.25) 



Also 7^],' 1 loop can be rewritten in the form 



(3), 1-loop 
/CT 



(m), 1— loop 
+7CT 



(5.26) 



where 7qt 1 l0 ° P formally corresponds to the topology (3) of Fig. 4, but now with the insertion 
only of the modified mass counterterm 



41 = m- 2£ (J) (£j 



ie 



The only nontrivial one-loop graphs with counterterms are therefore 7^' 1 loop and 7(jp ' 1 l0 ° P - 
The first one can be written as 



(1), 1— loop /-~<(X)i 1— loop 



UV 



+ 



1R 



+ c, 



(1), 1-loop ^(1), 1-loop 



(5.28) 



+ eCj J 



+ 0(e s 



where we have explicitly pointed out the IR divergent part. The individual orders of the e expansion 
are then 



^(1), l-loop _ 1 

°-i,uv - A 



.o(l), 1— loop 



(1), l-loop 







H(0, 0; z) + 2H{\, 0; z) - -vr 2 - 7^(0; z) 



+ m(H(0;z) + 2H(l;z))-iir/ 



-(7-3) 



(5.29) 
(5.30) 

(5.31) 



(1), l-loop 



H(2, 0; z) + -#(0, 0, 0; z) + #(1,0, 0; z) + 2#(1, 1, 0; z) 



-^y(H(0, 0;z) + 2H(1, 0;z)) 

1 77T 2 97T 2 

+ff 2 H(0; z) - ^H(0; z) - —H(l; z) 

+m \H(2; z) + iiT(0, 0; z) + ff(l, 0; z) + 2H{1, 1; z 

-^7(^(0; z)+ 2^(1; z)) 

i7T 3 7T 2 

T -C(3) + T 7 



ITT 



^2 



+ 4^ 



1 2 3 ?H 7 

H — 7 7 H 1- -■ 

8 ' 4 ' 48 4 



(5.32) 
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t-, (m), 1— loop , 

±<or 7ct we 



(to), 1-loop 
7CT 



C 



-4e ( f L 



a 



(to), 1— loop 
-1 



+ Q 



(to), 1— loop 




+ eC\ 



(to), 1— loop 



+ 0(e 2 ) 



(5.33) 



where 



(to), 1— loop 



(m), 1— loop 




(to), 1— loop 



1 

2^ 



#(0, 0; z) + i7riJ(0; z) - 2iT(-2; z) + 



7T 



-4H"(-1; z) + 2H(0; z) + 2ivr + 2 



1 

2p 



(5.34) 



(5.35) 



2H{-3; z) + AH (-2, z) + 4H(-2, -1; z) 

+4ff(l,-2;z)+4ff(-l,-2;js) 

-2#(0, 0; z) - 2H(-1, 0, 0; z) + ff(0, 0, 0; z) 

-2F(l,0,0;z) - 2iT(-2,0;z) 

-ivr(2#(0; z) + 2H(-1, 0; z) - H(0, 0; z) + 2ff(l, 0; z) + 2fT(-2; z)) 
-7 ( ff(0, 0; z) - 2H(-2; z) + iTrff (0; z) + y 

- z) + (-1; z) + 2tf (0; z) + 1) + 5C(3) - — 



+2 



4# (-1, -1; z) - 2H(-2; z) - 2H(-1, 0; z) + H(0, 0; z) 



2vr 2 



+(7 - l)(2ff(-l; z) - (0; z) - 1) - — + 1 



-i7r(2F(-l;z)-if(0;z)+7-l) 



(5.36) 



6 The two-loop graphs 

This section is devoted to the six 0(a 3 p 2 ) two-loop graphs depicted in Fig. 3. Let us denote their 
contributions to the re-scaled form factor 7 defined by (5.1) as 7W' 2 ~ 1o °p. Due to the symmetries 
of the graphs, we get the relations 



(2), 2-loop = (4), 2-loop _ (T), 2-loop 
(5), 2-loop = (6), 2-loop = (II), 2-loop^ 



(6.1) 
(6.2) 



We have therefore only four independent two-loop contributions to the form factor 7. For their 
determination we use the standard procedure based on the reduction to the scalar Master Integrals 
(MI) and on calculation of the latter by means of the differential equations technique. 



22 



6.1 Reduction to Master Integrals 

As a first step we define set of basic scalar integrals 

d d k 



B(m, . . . ,m) 



ic 



d d l 



n 



(6.3) 



(2vr) d (2vr) d AX A(M) n< 

where ni are integers and {Di(k,l)}J =1 is a set of seven independent propagator denominators, 
namely 



Di(M) 


= Z 2 -m 2 


£> 2 (M) 


= (Z + 9-) 2 




= (i - 9+) 2 




= A; 2 


£> 6 (M) 


= (A; + ^) 2 


£> 6 (M) 


= (^-<z + ) 2 


L> 7 (M) 


= (i-*) 2 " 



(6.4) 



m 



m 

2 



which at the same time form the basis for the seven independent scalar products build from the 
four-vectors k, I, q+ and q-. The momentum flow in B{n\, . . . ,n 7 ) corresponds to the auxiliary 
diagram shown in Fig. 5. The integrals B(ni, . . . , nj) depend besides the electron mass m only on 
one independent scalar variable Q 2 = (g+ + q~) 2 and as a consequence of the symmetry properties 
of the auxiliary diagram they satisfy the relation 



B{n 1 ,n3,n2,n^,n 6 ,n 5 ,n 7 ) = B(m, n 2 , n 3 , n 4 , n 5 , n 6 , n 7 ). 



(6.5) 



Note also, that some of B(n±, . . . ,n 7 ) are identically equal to zero. For instance whenever rij < 
for i = 1,2,7 (or i = 5,6,7) simultaneously, we can factor out a massless tadpole and therefore 
B(ni, . . . , nj) vanishes. 

Using the projection (2.8) and expressing the scalar products in the numerator of the inte- 
grands in terms of Di(k,l) we get schematically 

'a' 



7 



(i), 2— loop 



-4s 



,(0 



(m, . . . ,n 7 ;y)B(ni, . . . ,n 7 ) 



(6.6) 



711,. ..,717 



where cW(ni, . . . , n 7 ; y) are known coefficients which depend on the variable y = Q 2 /4m 2 and 
electron mass m. The explicit form of the reduction formulae (6.6) can be found in the Appendix 
B. The relations (6.2) can be explicitly verified for the right hand side of (6.6) using (6.5). 

There are altogether 172 integrals B(n\, ■ ■ ■ ,n 7 ) entering the sums in the reduction formulae 
(6.6), however, not all of them are independent. In addition to the symmetry property (6.5) there 
are also additional relations based on the integration by parts (IBP) [29, 30] and Lorentz invariance 
identities (LI) [31]. The former relations are consequences of the vanishing of the integral of the 
total divergence within DR. In our case we get eight relations schematically written as 

( v\ 



f d d k 


d d l i 


' 9 \ 


1 (2Tr) d 


(27T) d \ 


KW 1 J 



\<f-J 



n 



\ Di(k, l) ni 



0. 



(6.7) 
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Figure 5: The auxiliary diagram describing the momentum flow in the scalar integrals 
B(yi\, . . . ,nj). In this and the following Figures 6.-9., the full internal lines correspond to the 
massive scalar propagators with mass m while wiggly lines stay for massless ones. The full ex- 
ternal lines carry the momenta on the electron/positron mass shell and dashed external lines 
correspond to the momentum Q. 




Figure 6: The subset of MI with two propagators. 

The remaining LI relations express the invariance of the scalar integrals B{ri\, . . . , n-?) with respect 
to the Lorentz transformation of the external momenta 

'+/- - + f-lT- - 1-JT-) B ^ ■ - - >"7) = 0. (6-8) 

^dq +v ^dq +fl dq- v dq^^J 

The left hand sides of both (6.7) and (6.8) (when contracted with q+q^L) can be expressed in 
terms of linear combinations of B(ni, . . . , rij) with various rij and with y dependent coefficients. 
The explicit form of the resulting relations is postponed to Appendix C, here we give only the LI 
identity as an illustration: 

[(1 - 2y) (n 2 + n 5 ) - n 2 3~2 + + 4ym 2 n 2 2 + + 2yn 2 l~2 + - n 5 6~5 + + 2yn 5 4~5 + 
-(2h3, 5«6)]%,...,n 7 )=0, (6.9) 

where we have introduced the usual operators j ± (j = 1, . . . , 7) which act on B(n\, . . . , n-j) as 

j ± S(ni, ...,rij,... n 7 ) = B(m, ...,nj±l,... n 7 ). (6.10) 

The above linear relations (6.7) and (6.8) together with the symmetry property (6.5) can be 
used as an input for the Laporta-Remiddi algorithm [27, 28] which allows further reduction of 
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M (<0 (/) 



Figure 7: The subset of MI with three propagators. There are five types of different topologies, 
namely (a), (b), (c), (d) and (e, f). 




Figure 8: The subset of MI with four propagators. There are five types of different topologies, 
namely (a), (b), (c), (d, e, f) and (g, h). 

the number of the scalar integrals B(n\, . . . ,717) expressing them in terms of smaller number of 
MI. We have used the Maple implementation of the reduction procedure AIR [51] and found that 
all the necessary 172 integrals B{n\, . . . ,717) can be expressed in terms of 18 MI which belong to 
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Figure 9: The subset of MI with five propagators. There are two topologies within this subset, 
namely (a) and (b,c) 



12 different types of topologies, namely one two propagator MI (Fig. 6, one topology), six three 
propagator MI (Fig. 7, five topologies), eight four propagator MI (Fig. 8, five topologies) and 
three five propagator MI (Fig. 9, two topologies). In the above figures, the dots added to the 
internal propagator line mean that corresponding rij > 1 (number of dots is rtj — 1) while crossed 
line indicates that rii < 0, i.e. the i— th propagator is missing and the integral is a tensor one 
(number of crosses corresponds to \n>i\). The set of MI is summarized in Tab. 1. As a final 



Number of propagators 


MI 


2 


5(0,0,0,0,0,1,1)* 


3 


5(1,0,0,0,0,1,1), 5(0,0,1,1,0,0,1), 5(0,0,0,0,1,1,1)*, 
5(0,1,1,0,0,0,1)*, {5(0,2,0,0,0,2,1), 5(0,1,0,0,0,2,2)} 


4 


5(0,1,1,0,1,1,0)*, 5(0,1,0,1,0,1,1), 5(1,1,1,0,0,0,1)*, 

{5(1,1,0,0,0,1,1), 5(2,1,0,0,0,1,1)}, 

{£(-2,2, 1,1,0, 0,1), 5(0,1,1,2,0,0,1), 5(0,1,1,1,0,0,1)} 


5 


5(0,1, 1,0, 1,1, 1),{ 5(-l,l, 1,1, 0,1,1), 5(0,1,1,1,-1,1,1)} 



Table 1: The MI for different number of propagators ordered according to the figures 6-9. The 
simple products of one-loop integrals are denoted by star. The MI belonging to the same topology 
class are placed in curly brackets. 

result of the reduction procedure we get the individual two- loop contributions in the form of linear 
combinations of the MI with y dependent coefficients 
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According to the above formulae we shall need the MI up to (and including) the order 0(e). 
However, for the calculation of MI as described in the next subsection we need some of them to 
the order 0(e 2 ). 



6.2 Calculation of the Master Integrals 

Some of the MI are not genuine two-loop integrals and can be calculated simply as a product of 
one- loop integrals. This concerns four MI denoted by star in Tab. 1. The remaining fourteen MI 
can be obtained in a standard way as a solution of the appropriate closed system of ordinary linear 
differential equations [32, 33, 34, 35, 36]. The latter can be obtained by means of differentiating 
the integrals B(n\, . . . ,n 7 ) with respect to q±. On one hand, the result of such a differentiation 
contracted with external momenta can be related to the derivative of B(ni, . . . ,n 7 ) with respect 
to the variable y as 
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(here the second identity is a consequence of the LI). On the other hand, the right hand side of 
(6.15) can be expressed as a linear combination (with y dependent coefficients) of the integrals 
B(ni, . . . , rij) belonging either to the same topology class or to the topologies with less propagators. 
Finally we get (in terms of the operators (6.10)) 8 
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(6.16) 



Writing these equations for the MI we get at the right hand side along with the original MI also 
other B(ni, . . . , n 7 ys which can be subsequently expressed in terms of the MI to close the system. 
From the form of (6.16) it can be seen that the resulting set of differential equations has "triangular 
structure" in the sense that it connects the derivative of given MI either with MI with the same 
topology or with topologies with smaller number of propagators. 



It is easy to see that the difference of both possibie right hand sides of (6.16) just corresponds to the LI identity 



(6.9). 
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For further convenience we introduce dimensionless quantities b(ni, . . . ,n 7 ) denned as 

B(n u ...,n 7 ) = (iT(l + e){4Tr) £ - 2 ) 2 m 2 ( 4 "^ ni )b(n u ...,n 7 ) (6.17) 

and rewrite the differential equation for b(ni, . . . , n 7 ) in terms of the variable x (cf. (2.10)). The 
system (6.16) for b(n\, . . . , n 7 ) is then solved order by order in the e = 2 — d/2 expansion writing 

b(ni,...,n 7 ) = ^2 h{ni,. . . ,n 7 )e l (6.18) 
i>-2 

and expanding the right hand side of (6.16) to the given power of e. At each order, we solve 
the corresponding equations in the unphysical region < x < 1 where the MI are analytic. The 
solution is then fixed uniquely up to the integration constants which can be determined by the 
requirement of the absence of singularities in some appropriately chosen points of this analyticity 
region. For the calculation we have used the Mathematica package HPL [52, 53]. The results 
are summarized in Appendix D where also comparison of our independent calculations with those 
existing in the literature is given. 



6.3 Two-loop contributions 

Inserting the results of Appendix D in formulae (6.11-6.14) we obtain the final form of the two-loop 
contributions. Writing them in the form (cf. (4.8)) 
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+ ^ (367 (27 - 25) - 322. 2 + 1791) + 



-17T 
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24 
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17T 
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7T 
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(H(-2, z) + #(-1, 0, *) + #(1,0, *)) + — 





(6.31) 



The graph (1) has beside the UV also IR divergences, which can be identified using the general 
identities from Subsection 4.3. We get for j = — 1, —2 



7 



(1), 2-loop _ (1), 2-loop (1), 2-loop 



Tj,UV 



+ T 



i,iR 



(6.32) 
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where 



(1), 2-loop 
7-2, UV 

(1), 2-loop 
7-2, IR 

(1), 2-loop 
UV 



(1), 2-loop 
7-1 IR 



_3_ 
16 



-[/3 + -)(H(0;z)+m) 



5 



1 + 



1^ 



H(0, 0, z) + 2H(1,0, z) + ivr (H(0, z) + 2H(1, z)) - -vr 2 



7(^(0,z)+i7T) 



3 35 

+ -7 

8 32 



2H{-3, z) + (-2, 0, z) - -H(0, 0, 0, z) 



ff(-2,z)-|ff(0,0,*)-^ 



1 



1-77 /9+o +"5- 1 + 1* 



7T" 



6 



1 



(H(0,z)+m) 



(6.33) 
(6.34) 

(6.35) 



(6.36) 



From the above formulae the general relations (4.12), which are valid for each graph separately, 
can be easily verified. 



6.4 Two-loop counterterm contribution 

For the construction of the two-loop counterterm (3.6) we need further the local parts of the 
0(e~ 1 ) UV divergences of the two-loop graphs (cf. (4.11)). Using the results of the previous 
(sub)sections we get 



(1), 2-loop\ 

7-i, uv ; 
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l— loop (n), l— loop 
7o 7li 
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" ~32 

/~<{m), 1— loop (m), 1— loop 

Co 7-i 
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(6.37) 



where 7^' 1 loop are the coefficients of the UV divergent parts of the one-loop (sub)graphs (see 
(4.3)). According to general formula (4.15) we get for the two-loop counterterm contribution 



tree 

7cr 



m 



ot\ 



e + 4(7 



(1), 2-loop 
' -2, UV 



+ 2 7 . 



(r), 2-loop (3), 2-loop 



(n), 2-loop \ 



^(ci 1 i ) ; u - v loop x + ^ 



x 6 - 7 :T loop ^ + 
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(3), 2-loop\ 
h I7_i, UV J 



+ 7-2" " ^ + 2r- 2 " ' ) 

) 



+ 2 , (n), 2-ioop^ 



and for our choice of the renormalization scheme 



tree 
7CT 
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m) Kir) 
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^ + 8^ " e U X ~8 7 + l8" 



where 
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In 



in. 



HI 

ni- 



+ 0(e). 



(6.38) 

(6.39) 
(6.40) 
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To get the numerical values of the two- loop radiative corrections, we need to know the NLO coun- 
terterm coupling £ or its scale dependent renormalized value £ r (/i). In principle it can be obtained 
similarly as x r (^) ( cr> - [6]) by means of matching of the (complete) NLO chiral expansion of the 
amplitude with the sum of the same types of graphs as we calculated but now with appropriate 
model of nonlocal off-shell pion transition form factor F 7r o 7 * 7 * in place of the local 7r°77 vertex 
(3.2). This is however beyond the scope of our paper. Instead we make a simple estimate of the 
value of £ r (/i) using its running with the renormalization scale. From (6.40) with central value of 
X = —16.8 (see (5.11)) we get 

AC = If(lGeV) - f(0.5GeV)| = 5.5 (6.41) 
We therefore roughly estimate 

f(770MeV) = 0±5.5 (6.42) 

and get finally 

|= -32.3 + 3.7(x + 16.8) ±5.5. (6.43) 



7 Soft photon bremsstrahlung 

Within the soft photon approximation, the amplitude of the process tt° — > e + e~7 factorizes 

'(g_- £ *(M)) (q+-e*(k,\)y 



(7.1) 



(q- ■ k) (q+ ■ fe) 

where e(k, A) is the polarization vector of the emitted photon with soft momentum k and helicity 
A. Taking the square of the modulus, summing over helicities and integrating over the soft photon 
region |k| < u defined by the experimental energy cut 



we get 



where 



r 7r o„ i , e + e - 7 — Ibs^-h-o 



7r u — >e~T~e~ 



Ibs 



d 3 k 



k , <w (2vr)32|k| 



m 



m 



(7.2) 

(7.3) 
(7.4) 



.(<Z- " k)(q+ ■ k) ( g _-fc)2 (q + -k)\ 
The latter integral is IR divergent and can be regularized using DR. Let us divide the resulting 
Ibs into the diagonal part 



.diag_ _ ll -2e(]L\ 2e (°\ [ d3 2 " k 

BS ~ " W U J w<u (27r)i-*2|k 



m 



2+2e 



+ 



m 



2+2e 



{q.-kf (q + -ky 

d 3 " 2e k 2m 2£ (q + -q_) 
.n J 7, k | <w (27r)i-2^2|k| (q--k)(q+-ky 

In the rest system of the decaying pion both integrals are elementary and we get 



and non-diagonal part 

r non— diaa 



I 



BS 



mJ 



(7.5) 



(7.6) 
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L BS 
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2e(V\ 2£ ( a \ ( 47r 



m 



ttJ T(l-e) Iw 



in 



2= 



1 1 

e ~ 



H(0;z) -ln4 + 0(e) 
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and 



/ in \ 

---#(0;z)-ln4 + 2m(-) ~7 + 0(e) 



riion- diag 
l BS 



\mJ \irJ 2T(1 - e) \oj 
H(0;z) 



7T 



e 



2 In 2F(0; z) - 2H(1, 0; z) - #(0, 0; z) + 0(e) 

3 



m 



^^1-2H(1,0;z)-H(0,0;z)-y 



CO 



+H(0;z) 21n — -7- In 4 +0(e 



(7.7) 



(7.8) 



Note that the IR divergent part of I B l g g coincides up to a sign with that of the 2dT,(p)/dfl (see 
(5.21)). The latter factor is necessary for the renormalization of the external fermion lines. Namely, 
on the level of the decay width 



Z 2 e T^ e+e - = ( 1 + 2^g^ + 0{a 2 ) ) I>_ 



(7.9) 



As we have mentioned in Subsection 5.2, within our renormalization scheme 

of \mJ \ttJ e 

Therefore up to the assumed accuracy we can effectively make the following replacement 



(7.10) 



(7.11) 



provided we replace I^g S with its finite part. 

Taking this modification into account we can finally write the bremsstrahlung integral Ibs m 
the form (4.25) 



Ibs = V 



mJ V7r 



— + I + 0(e) 



with 



/ 777, \ 2 

/o = -^ ( 0;z) + ln(-) -7 



(7.12) 



(7.13) 



(7.14) 



The relation (4.27) is now manifest. This completes the list of all the necessary ingredients for the 
final calculation of the radiative correction 5(x^ t ). 
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8 The two-loop radiative correction 



8.1 The exact two-loop result 

Putting the results of the previous sections together and using the general formula (4.31) we can 
finally write for the complete QED two-loop correction (2.13) schematically 

S ( x ff) = j + Re (_L ( 2 £ + A« + 2A< r > + A< 3 > + 2 A™)) 

= § bs + + + ^r) + 6 (3) + § {U) (gl) 



In this formula 7^0 = 7o ~ loop + \ and the variable x c ^ is connected to the maximal energy 
of the soft photon included into the inclusive 7r° — > e + e~7(|k| < 00) decay rate by (7.2). In the 
above expression we have explicitly separated the contribution coming from the bremsstrahlung, 
the RG invariant NLO coupling £ and the individual graphs: 

A« = 2 7 W' 2 " loop + 2xC? } ' ^ loop + 3Cf >• 1 ~ loop - 27 1 1 " loop ci 1 1 ) ; I ^ loop 

A< r > = 2 7 S r) ' 2 - loop + 2x 67o 1 - loop - ^-^p 

= 27 (D, 2 -ioo P+ 3^__0 7ol _ loop _l 7il _ loop 

A( 3 ) = 2 7 J 3) ' 2 " loop - 2X670 1 " 10015 + 2(x 6 + x 7 )C^ 1 " loop + ^l~ loop + ^c{ m) ' ^ loop 

_ o„,( 3 )' 2-loop ^ / *M / 1-loop . n {m), l-loop\ 1 1-loop 3 „(m), 1-loop 

" 270 2 V 7 V V 7 ° + °° J + 2 71 + 2 1 

A< n > = 2 7 ^ n) ' 2 " lo ° P + 2x 8 7 ^ loop + ^ 7l - Ioop 

3 

(II), 2— loop 2_ i_loop . ^ i_loop /o \ 

= 27o - 377o +371 • ( 8 - 2 ) 

Here we have inserted for X{ the particular values corresponding to our choice of the renormaliza- 
tion scheme (cf. Section 5.2) and 7q*' 1 ' 2 loop , Cq 1 loop , c[^' 1 loop , C_j jp loop and 7^~ loop are 
explicitly given in Subsections 6.3, 5.3 and 5.1 respectively. 

The numerical results for various contributions are summarized in the second column of Tab. 
2. Here we use for the constants x an d £ the values (5.11) and (6.43) and the following numerical 
entries: M w a = 135MeV, m = 0.51MeV and a = 1/137. For the sum of all the contributions we 
get finally 

$G) + + «5( r ) + j(3) + j(n) = (_ .8 ± 0.2) %, (8.3) 

where as the only source of errors we take the uncertainties of x an d £■ We observe considerable 
cancellation between the large contributions 6^ and 5^ which makes the role of the relatively 
smaller contributions of the other graphs numerically important. 
Adding the bremsstrahlung we have 

6(xff) = (4.71n(l - x^ ut ) + 8.3 ± 0.2) % (8.4) 

and for x c p = 0.95 which is the cut used by KTeV we get 

S{x^ = 0.95) = (-5.8 ± 0.2) %. (8.5) 
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8.2 Large-logarithm approximations to the exact result 

The exact two- loop expressions for 5(x c ] ^ t ) and for are rather long but they can be approxi- 
mated with a very good accuracy performing the large-logarithm (LL) expansion in terms of 



L 



Inz 



In 



1-/3 
1 + /3 



I 7T U 



(8.6) 



and taking into account only the leading terms (up to the order 0{L For the various 

contributions we get 9 
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9 In these and following formulae we omit also all the terms of the order 0(1 — f3). 
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( '[7o\ 


exact 


T T 4-" 1 

LL, rational 


T T 1 "1 

LL, polynomial 


LL, polynomial, KG 


X -> oo 


iff 1 ) 


-0.36 


—0.36 





—0.92 







18.4 


18.4 


18.8 


19.5 


—7.29 


x(D 
( ' 


—22.2 


—22.2 


—24.6 


—22.8 





A(3) 


n no 


n no 


Q A Q 


Loo 


n 
U 




2.45 


2.45 


1.17 


2.41 





5(x c ^ t = 0.95) 


-5.84 


-5.84 


-6.29 


-5.28 


-12.4 



Table 2: Numerical values of the contributions of the individual graphs to the total two-loop 
correction 5 within various approximations described in the main text. 
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5.13) 



Inserting the above expressions into (8.1) we get an approximation of 5(x ( j^ t ) in terms of the 
rational function of the variable L. It would be tempting to expand further the factor i R 
(8.1) and approximate the whole 5(x c ] ^ t ) as a second order polynomial of L with the result 



^LL polynomial ^cut 1 

or numerically 
and 



^ + ^L + 2(L-l)ln(l 



13. 
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S LL polynomial ^cut) = (47^ 



~ C f) ) - -x - —1 - + - 
D ' 3 A 2 ' 36 9 



+ 7.9) 



§ LL polynomial ^cut = Q 95) = _g ; 



.14) 



5.15) 



5.16) 



We can also proceed apparently more carefully and include part of the large RG logarithms into 
the expansion writing In (^ 2 /m 2 ) = — L + In (/i 2 /M 2 ); such a modification makes a difference in 
the 0(L ) terms of the expansion. We get in this case 10 



s-LL polynomial / cut \ 
RG \ X D ) 



L 8 r ntr f1 „ t . 1 13 19 43 

h -L + 2(L - 1) ln(l - x*?r) y 7 vr 2 H 

12 9 v ; v D 7 3 A 2 ' 36 9 



3 In the final expression we expressed again x r (M) an d ln(/x 2 /M^ ) in terms of x and £■ 
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(4.71n(l 



x^) + 8.9) 



5.17) 



which gives 

^LL polynomial (3 ,cut = Q %) = (81g) 

In both cases the reason for not very good agreement with the exact result is that the expansion 
of I/7L0 does not converge well 11 and therefore it is much safer to approximate 5(x ( j^ t ) in terms 
of rational function of L. We have compared all the possibilities (LL rational, LL polynomial and 
LL polynomial with RG logs included) in the third, fourth and fifth column of Tab. 2 respectively. 
The rational LL approximation is in excellent agreement with the exact result. 

The LL approximation was originally calculated by Dorokhov et al. in [21] including only 
the graphs (l)-(4) in the Fig. 3 and taking into account only selected regions of the two-loop 
integration space. As we know from (8.12), the omitted graphs contribute only in the next to 
leading order in the LL expansion, however due to the accidental cancellation of the leading order 
terms they are in fact numerically important. The result of [21] for the virtual and soft photon 
corrections is 



Dorokhov 
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24 



L 2 + 2(L-l)ln (l 



™cuf 
X D 



7T- 



+ 2 



(8.19) 



(4.71n(l 



*D Ut ) + 6.7) 



which numerically gives for x^ = 0.95 a correction with significantly larger absolute value in 
comparison with (8.5), 

C^hof 7 « = 0.95) = -13.3%. (8.20) 

The formula (8.19) could be compared with our polynomial LL approximations after subtracting 
the contributions of the graphs (5) and (6) in Fig. 3 (i.e. those with vacuum polarization insertion 
into the internal photon lines). In the two variants described above we get however a result 
substantially different from S^^''j g f°£ t 7 (x™ t ), namely 



rLL polynomial / cut 
'(l)-(4) 



and 
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and numerical values for x^ 



0.95 

LL pol 

'(l)-(4) 



rLL polynomial / cut 
On-* ia\ [X D 



rLL polynomial / cut 
or* /i\ I a\ \X£) 



0.95) 
0.95) 



(8.21) 

(8.22) 

(8.23) 
(8.24) 



RG, (l)-(4) 

The reason of this discrepancy is difficult to trace out because a completely different framework 
has been used for the calculations in [21] and we therefore left this problem open for further study. 



"While 1/jlo = 0.024 + 0.044i, the approximation up to the order 0(L" 6 ) gives [1/7lo] (5) = 0.044 + 0.037L 
Including part of the large RG logarithms into the expansion as described in the main text we get [I/7L0]' 5 ' = 
0.047 + 0.053L 
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8.3 Point-like n°e + e vertex approximation 

Let us briefly comment on another type of approximation which is connected with the model calcu- 
lation of the QED radiative corrections by Bergstrom [12]. His result including virtual corrections 
and soft photon bremsstrahlung reads 



virt.+soft 7/ cut^ ( a 



"Bergstrom \ D ) ~ 



2(L - 1) In (1 - + y - 1 + 0(1 - xff) 



(8.25) 



where the 0(1 — x c p l ) terms, which we do not write explicitly, stem from the fact that the real 
photon radiation was calculated exactly i.e. beyond the soft photon approximation. Taking also 
these terms into account we get numerically 

CrgsS 7 « = 0-95) = -13.8%, (8.26) 

which is remarkably close to (8.20). However, the calculations [12] used a completely different 
approximation from the LL one that was used in [21]. 

In the paper [12] the approximation was based on the substitution for the nonlocal one-loop 
7r°e + e _ (sub)graph with a local effective vertex of the form 

£ c ff = ifiU^evr . (8.27) 

In appropriate renormalization scheme this vertex is essentially equivalent 12 to the vertex 

Cos = -f^7 M 75e<V° = -^r )' Xcsn^ed^ , (8.28) 

which takes properly into account the GB nature of pion and has the same structure as the 
counterterm (3.3). The Bergstrom's calculation can be therefore qualitatively understood as the 
leading order term in the formal large x expansion of the full twodoop result 13 (i.e. corresponding 
to the assumption of small nonlocal part of onedoop 7r°e + e~ (sub)graph in reference to its local 
part represented by x). Performing further the LL expansion of this leading term we get for the 
Bergstrom-like approximation of our result 1 



7T 



2(L - 1) In (1 - x c D ut ) + ^ + ^(l- 7 ) + {xr\L- x ) 



129) 



For illustration purposes we add the corresponding numerical values as the sixth column of Tab. 
2. However for physical value of x such an expansion does not converge (note that |7o~ loop /x| ~ 2) 
and therefore the Bergstrom's result can not be taken as a serious approximation of the full two- 
loop d(x^ t ). 



9 The phenomenological applications of the results: first look 

In this section we discuss several issues connected with the phenomenological aspects of the results 
obtained above. Up to now we have fixed or estimated the free parameters x r (^) an d £ r (/-t) of 

12 We ignore here the axial anomaly which does not contribute to the relevant (e + e~ |9 M e7 M 75e|0) matrix element. 

13 Note that in the limit % — > 00 only the contribution of the one-loop graph (1) of Fig. 4 and soft bremsstrahlung 
are effectively taken into account. 

14 The difference between (8.25) and (8.29) corresponds to different regularization of the IR divergences and 
different renormalization scheme. 
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CLEO bound 
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CLEO+OPE 


QCDsr 


LMD±V 


QM 


NxQM 


VM 
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— 17.8 


— 16.5 ± 0.3 


— 16.3 ± 0.1 


—16.5 


—18.0 ± 0.5 


— 16.7 ± 0.5 


— 19.1 


x r 


1.3 


2.6 ±0.3 


2.8 ±0.1 


2.5 


1.1 ±0.5 


2.4 ±0.5 


-0.05 


6[%) 


-5.68 


-5.89 


-5.92 


-5.88 


-5.65 


-5.85 


-5.50 



Table 3: Illustration of the sensitivity of the central value of the two loop QED correction 5 on 
the various values of x described in the main text. The renormalized \ r is taken at [i = M p . 




-4 -2 2 4 



Figure 10: The dependence of 5(x ( j^ t = 0.95) on x r (M p ). The filled band corresponds to the 
variation of 5(x'j^ t = 0.95) with ^ r {M p ) inside its error bar. 

the effective Lagrangian and concentrated on the corresponding prediction of the two-loop QED 
corrections. In such a way obtained 5{x c ^ l t ) represents a part of the theoretical tools necessary for 
a precise extraction of the 5 no - rad ( 7r _ e + e~) from experimental data. Therefore a consistency 
check of our final result (which uses one specific value of x r (/ i )) with that using other estimates 
of x r (^) which are available in the literature should be desirable. Such a check is the topic of the 
first subsection. 

We can also, however, reverse the point of view and investigate the sensitivity of the theoretical 
prediction for the branching ratio B(ir° — > e + e _ (7),X£> > 0.95) on the free parameter x r (A*) an d 
try to extract the information on its actual value from the experimental data. Though we have 
not all the necessary ingredients at hand, we can make a preliminary analysis of this issue. This 
is done in the second subsection. 

Last but not least, our result is closely related to the calculation of two-loop QED corrections 
to the class of processes of the type P — > l + l~ where P = tt°, i], Kl and I = e, /j,; this relation is 
briefly discussed in the last subsection. 

9.1 Note on the dependence on % 

In the above numerical calculations we have fixed the value of the constant x according to the large 
Nc inspired LMD estimate (3.4) of the effective coupling x r (p) entering the Lagrangian (3.3). In 
the literature there exist, however, further model dependent estimates of this constant based on 
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various models or phenomenological parameterizations of the pion transition form factor F w o~*~* 
(see [11] for comprehensive review). In Tab. 3 we summarize various values of % and x r (M p ) 
which we take over 15 from [11]. The first three columns (denoted as CLEO bound, CLEO+OPE 
and QCDsr) correspond to various treatments of the parametrization of F n o~*~* (t, t) using the 
CLEO [10] data (see [11] for details), next column (LMD+V) is the improvement of the large 
Nc estimate mentioned above including two 1 multiplets 16 [6]. The column QM is based on 
the constituent quark model while the last two columns (N%QM and VM) on two variants of 
the nonlocal chiral quark model [54] and [55] respectively. In the the last row we illustrate the 
sensitivity of our result for the central value of 5(x c ^ t = 0.95) on %. All but the last resulting 
central values of 5{x c ^ t = 0.95) are compatible with LMD result (8.5) within the estimated error 
bar. The dependence of 5{x c ^; t = 0.95) on x r (M p ) in a wider range is plotted in Fig. 10 where 
also the variation with £ r (Mp) inside its estimated error bar is illustrated by the filled band. ' 

9.2 Note on the phenomenological determination of x r (M p ) from ir° — > e + e~ 
decay 

Let us stress that the above analytical result for 5(x ( j^ t ) represents only a part of the problem 
of the complete QED radiative corrections to the process under consideration, and that the 
realistic analysis of the experimental data requires several additional pieces of information. The 
first one corresponds to the issue of the hard photon bremsstrahlung for which the soft photon 
approximation is not an adequate framework and for which more appropriate calculations have 
to be done. A closely related issue is the applicability of soft photon approximation for the 
KTeV choice of the cut x c ^ 1 = 0.95. Another missing information is connected with the Dalitz 
decay which yields the same final state as the real photon bremsstrahlung in the 7r° — > e + e~ 
decay. Though this process is dominated by low xd and is therefore suppressed for xp > 0.95, its 
integrated contribution is known to grow rapidly with decreasing x^. 

These additional issues have been addressed in the present context already in the paper [12], 
(see also [21]) and in this form they have been used for the analysis of the experimental data by 
the KTeV collaboration [8]. However, the analysis performed in [12] might be incomplete. As 
far as the hard photon bremsstrahlung is concerned, the point-like 7r°e + e _ vertex approximation 
described in Subsection 8.3 has been used. This approach is, however, well justified only for the 
soft photon region where the details of the off-shell 7r°e + e _ vertex are inessential. 

Also, the corrections due to the Dalitz decay calculated in [12] (based on the radiative cor- 
rections to 7T° — > e + e~7(7) obtained in [56]) and used in [8] should be taken with some caution. 
As it has been shown recently [57], the one-photon irreducible (ljIR) contributions which were 
omitted in [56] are in fact important already for xq > 0.6 where they give a negative contribution 
5 llIR (xo) < —1% of the leading order differential decay rate (dr Dalltz /dx£>) L °. 

The importance of the detailed knowledge of dr Dalltz /dx£) is twofold. On one hand, because 
the Dalitz decay has been used by KTeV-E799-II as a normalization and because the measured 

15 In [11], the values of -4(0) = x + §7 — § are presented. 
16 This ansatz is denoted gVMD in [11] 

17 Note that for fixed \ the dependence of 8{x c q 1 = 0.95) on £ is trivial (i.e. linear as can be seen from (8.1)). 
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quantity was the ratio 

r(vr° -> e+e~,x D > 0.95) 



'1.685 ± 0.064 ± 0.027) x 10 -4 , (9.1) 



r(vr° -> e+e-7,x D > 0.232) 

it is necessary to extrapolate the Dalitz branching ratio to the full range of xo. As this ex- 
trapolation is concerned, the missing corrections to dr Dalltz /d:r£) are in fact inessential, since the 
contribution of l^IR corrections integrated over the full phase space can be shown to be one order 
of magnitude smaller then the experimental error of the r(7r° — > e + e _ 7)/T(7r° — > 77) branching 
ratio (see [57] for details). On the other hand, dr Dalltz /dxD has been used in order to subtract the 
Dalitz decay background in the region xd > 0.95. Here the effect of missing contributions might 
be more important. 

Let us give here only a preliminary illustration of the interrelation of our partial result of 
the QED radiative corrections and the precise branching ratio (1.1) obtained by KTeV-E799-II. 
Taking the above theoretical uncertainties into account we can write our prediction for the KTeV 
measured branching ratio as 

B(tt° -> e + e~(rf),x D > 0.95) = £(71-° -> 77) 



r(7r° —7- 77) 



1 + 5(0.95) + A B5 (0.95) + A l7/i? (0.95)), (9.2) 



where the only experimental input is the precise branching ratio B(ir° -> 77) = (98.823 ±0.034)%. 
In the above formula 

A BS (x^) ee 5 B x s act (x^) - 5g%(x^) (9.3) 
is the difference between the soft photon and exact bremsstrahlung correction and 

1 r 1 /jpDalitz \ LO 

= r^o^-) 4, <^ (9.4, 

corresponds to the unsubtracted fraction of the Dalitz decay background discussed above. Without 
detailed knowledge of the exact bremsstrahlung we can only roughly estimate the error A BS (x c ] ^ t ) 
of the soft photon approximation with help of the point-like ir°e + e~ vertex approximation used in 
[12]. We get 

A BS (x^) = -2 (?) (L - 1)(1 - * c D ut ) + O (1 - x^f) (9-5) 



. 7T 

„cut 



which gives for x™ = 0.95 a reasonable difference 

A B5 (0.95) ps -0.25%. (9.6) 

As far as the unsubtracted fraction of the Dalitz decay background is concerned, we can use the 
explicite formulae for the corrections taken from [57] and arrive at 18 

A^(0. 9 5)^- [rLO( ^; ;; MeV1 (9.7) 

which gives for x r (M p ) = 2.2 

A^ IR \ x r {Mp)=2 . 2 = -0.35%. (9.8) 



s Here we neglect a very weak dependence on the constant x r (M p ) in the numerator. 
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Figure 11: The dependence of predicted B(ir° — > e + e~(7), xd > 0.95) on x r (M p ). The filled band 
corresponds to the variation of the additional contributions A BS + A^ /IR from zero to its maximal 
estimated value described in the main text (the dashed line here represents A BS + A 1,yIR = 0). The 
dash-dotted line shows the leading order value. The horizontal band corresponds to the KTeV- 
E799-II measurement. The dashed and dotted vertical lines delineate the large Nc inspired LMD 
estimate [6] and the region compatible with experimental value respectively. 



Note that both these additional contributions are larger then the variation of 5(0.95) with £ r (M p ) 
inside its error bar (which yields A^(5(0.95) < 0.15%). 

In the Fig. 11 we have plotted the right hand side of (9.2) as a function of x r (M p ) both with 
and without the estimated additional contributions A and A llIR together with the leading order 
branching ratio against the experimental value (1.1). The vertical bands correspond to the large 
Nc inspired LMD estimate x r (M p ) = 2.2 ± 0.9 and to the range of values compatible with (1.1). 
These bands do not overlap, the preferred region of x'i^p) is shifted towards higher values 19 , 
namely x r {M p ) « 4.5 ± 1.1 which corresponds to 5(0.95, x r (Af p ) = 4.5) = -6.2%. 

9.3 Generalization to the P — > l + l~ decays 

From the point of view of X-PT the coupling constant x r {^) is universal 20 , because it enters the 
SU(3) generalization of the counterterm Lagrangian C ct given in Appendix A and is therefore 
connected with other processes of the type P — > l + P where P = tt°, r), Kl and I = e, /i. These 
decays are governed at the leading order by the following effective long-distance Lagrangian 



l-f-fldpp 

(9.9) 



19 Quite interestingly, even higher values of y r (M p ) have been obtained from similar analysis of the related decays 
Kl — > ft + [f~ , namely Y r (M p ) = 8.07 ± 0.20 or 5.84 ± 0.20 [58], and the bigger solution obtained from 77 — > [i + pi~ , 
namely y r (M p ) = 8.0 ± 0.9 [59]. Here the solution for \ r (M p ) shows two-fold ambiguity which is present also in the 
7T° — > e + e~ case, however we have fixed it keeping only the solution closer to the large Nc prediction. 

20 The analogous statement is true also for ^(a*)- 
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Figure 12: On he left: the x^ 1 dependence of the (partial) two- loop QED corrections 5(x c ] ^ t ) to 
rj l+r decays for x P (M p ) = 2.2 (solid line), x r (M p ) = (doted line) and x r (M p ) = 5.5 (dashed 
line). On the right: the x r (M p ) dependence of 5{x c ^ t = 0.95). The filled band corresponds to 
variation of £ r (M p ) inside its error bar. 
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Figure 13: The same plots as in Fig. 12 for the Kl — > l + l decays. 
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P — V / + / — 


TV — r e e 


Tj — y 6^~C 


V MA 4 


J\p \l jl 


Kp —?• e + e 


5(w = 3.37MeV)[%] 


-5.8 


-16.3 


-3.9 


-3.0 


-15.1 


S(x c ^ t = 0.95) [%] 


-5.8 


-7.8 


-2.5 


-1.8 


-7.3 



Table 4: Illustration of the typical values of the (partial) two-loop QED corrections for processes 
P —7- l + l~ with the same cutoff uj for the soft photon energy and with the same cutoff on i/j. 

where Cp is a normalization factor and Xpif 1 ) 1S t ne effective long distance coupling which can be 
split into its universal xPT and specific short-distance components 

xW = x r (^) + x s P d - (9.10) 

The latter is nontrivial only for the P = Kp case where it is well known [60] , [61] (see also [59] 
and [62]) 

x f L = -1.82 ±0.04. (9.11) 

The La grangian pi+i- is up to the normalization and shift in x r {f Ji ) identical with that we 
have used for the calculation of the 0(a 3 p 2 ) corrections to 7r° — > e + e~ decay which means that 
our general result can be almost straightforwardly used for other P — > l + l~ processes by means of 
the substitution M w o — > Mp , m — >■ mi and x r ~^ X r p- More precisely, in such a way we obtain the 
contribution of the graphs depicted in Fig. 3 where all the fermion lines correspond to the final 
state lepton flavour, that means that for the decays P —> we miss the vacuum polarization 

insertion graphs with electrons inside the loop. The results of such a generalization are illustrated 
in Tab. 4 where the corresponding corrections are listed (either for the same cutoff uj = 3.37MeV 
on the energy E 7 of the soft photon or at the same cutoff = 0.95 onxp = mf +l _ /Mp) and 
compared with the case of 7r° — ?• e + e~ decay. To obtain these numbers we have fixed x r (M p ) = 2.2. 
As another illustration of the typical values of the (partial in the sense of the missing graphs) two- 
loop QED corrections we have plotted their xf^ and x r (M p ) dependence in Figs. 12 and 13. 

10 Summary and conclusion 

In this paper we have calculated the two-loop 0(a 3 p 2 ) QED radiative corrections to the rare 
decay 7r° — > e + e~ including all the relevant graphs. As a result we have obtained exact analytical 
expression which takes into account the virtual photon contributions without any approximation. 
The IR divergences has been treated including real soft photon bremsstrahlung. The latter has 
been calculated within the soft photon approximation and added to the inclusive decay rate 7r° — >• 
e + e~("f). We have worked in the framework of the xPT with dynamical leptons and photons and 
parameterized the missing information on the details of the pion transition form factor in terms 
of two a priori unknown RG invariant couplings x and £ which also incorporated the large RG 
logarithms. The numerical value of the first of these couplings has been obtained using the analysis 
of [6] based on the large Nq matching and LMD ansatz, while the value of the second one (on the 
value of which our result proved to be much less sensitive) has been estimated from the running 
of the corresponding renormalized coupling (, r (fJ>) with the RG scale \i. 
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The main motivation of our work was to test the validity of two approximative calculations of 
the QED radiative corrections already existing in the literature, namely the Bergstrom's point-like 
7r°e + e _ (sub)graph approximation [12], which has been used for the analysis of the experimental 
data by the KTeV collaboration, and the more sophisticated large-log approximation [21] which 
claimed to confirm the applicability of the previous one. 

We have identified the Bergstrom's calculation as a leading term in the large x expansion of 
the exact two- loop result. We have found that this expansion did not converge and therefore could 
not be trusted without reservation. The numerical discrepancy between the exact result and the 
Bergstrom's one seems to confirm this conclusion. 

Next we have discussed several variants of the LL approximations derived from the exact 
two-loop result. The two polynomial ones differ from the exact result by roughly twice the error 
estimated from the uncertainty of the couplings x an d £ while the rational one gives an excellent 
agreement. However we did not manage neither to confirm the LL calculation of [21] nor to 
reveal the reason of its large difference from our result. The approximation [21] was obtained in 
completely different framework and we therefore left the final resolution open to further studies. 

Our final result numerically reads (see the definition (2.13)) 

5{x c ^) = (4.71n(l - x<b ut ) + 8.3 ± 0.2) % (10.1) 

and for the cut xffi = 0.95 chosen by KTeV 

5(x c ^ = 0.95) = (-5.8 ± 0.2) %. (10.2) 

Here the error stems from the uncertainty of x and £. Our result significantly differs from the 
previous approximative calculations [12, 21] 

«KKS W = °- 95 ) = " 13 - 8 %' C'o^f 7 (*2? = 0-95) = -13.3%. (10.3) 

and the change is in the right direction towards the agreement of the experimental data with the 
SM prediction. 

Let us note, that for the realistic analysis of the experimental data it is necessary to discuss 
carefully two further topics, which we have only partially included in this work, namely the real final 
state radiation beyond the soft photon approximation and the incorporation of the Dalitz decay 
contribution. The former is important because it is well known that the soft photon approximation 
is not much reliable except of very limited region in the phase space due to the small electron 
mass. The latter process, which has been recently revisited in [57], is known to yield a non- 
negligible background of roughly 2 — 3% [12] near the cut x^ = 0.95 chosen by KTeV. The 
more detailed analysis of these issues is still in progress. In this work we have performed only 
preliminary simplified analysis and found that the SM prediction for the branching ratio B(7r° — > 
e + e~(^),X£> > 0.95) might be reconciled with the experimental value for x r (M p ) ~ 4.5 ± 1.1, 
which is however off the predictions for x r (M p ) based on the phenomenological models of the pion 
transition form factor. 
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A The xPT Lagrangian with dynamical photons and leptons 

In this appendix we summarize the relevant parts of the xPT Lagrangian with dynamical photons 
and leptons. Following the notation of [25] and using the SU(2)l x SU(2)r variant of the theory 
(see also [63] ) we need the following terms of the complete Lagrangian 

£ = £wzw + 4- 7 + 4- 7 + 4? (A.l) 

where 



wzw 



N c 
32vr 2 



T y,vpcr 



C 



(2) 
e— 7 



\f^ v + e(i 7 ^ - m)e + \ (d ■ A)' 

1 /a s 



4- 7 = (^) xeei^D^e + (^) x 7 mee + - (^) xgF^F^ 
4? = (^) 2 e7 M 7 5 e [ X i{Q 2 {D^UU + - D^U+U) + X 2(U + QD fl UQ - UQD^U+Q)] 
and where 



(A.2) 



u 


i / 7T , 

= eXP ^( v V2^, 




Q = diag 


(I 


1 

~3 


DU 


= dU - irU + iUl 


s = U A 


-dU 






r 


= r--(r), T-- 


-i-5<0 








V 


= d^Vy - d u v^ - i 











The coupling x is connected with xi.2 according to 

1 



X 



(xi + Xi) 



(A.3) 



(A.4) 



7 V 



B Reduction to scalar integrals 

In this appendix we give the reduction of the individual two-loop graphs to 172 scalar integrals. 

2-loop = o ( 8?r ) 4 (&\ ^-4e 
m 2 \ttJ 

x [-64m 6 y 3 5(0, 1, 1, 1, 1, 1, 1) + 32m d y 2 B(0, 1, 1, 1, 1, 1, 1) 
+32m A y 2 B(0, 0, 1, 1, 1, 1, 1) + 32m A y 2 B(0, 1,0, 1, 1, 1, 1) 
-16m 4 yB(0, 0, 1, 1, 1, 1, 1) - 16m 4 yB(0, 1, 0, 1, 1, 1, 1) 
-4m 2 yB{0, -1, 1, 1, 1, 1, 1) + 8m 2 yB(0, 0, 0, 1, 1, 1, 1) 
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-4m 2 y5(0, 0, 1, 1, 1, 0, 1) - 4m 2 y5(0, 1, -1, 1, 1, 1, 1) 
-4m 2 y5(0, 1, 0, 1, 0, 1, 1) + 8m 2 yB(0, 1, 1, 1, 0, 0, 1) 
-4m 2 y5(0, 1, 1, 1, 0, 1, 0) - 4m 2 yB(0, 1, 1, 1, 1, 0, 0) 
+2m 2 B(0, -1, 1, 1, 1, 1, 1) - 4m 2 5(0, 0, 0, 1, 1, 1, 1) 
+2m 2 B(0, 1, -1, 1, 1, 1, 1) + 5(0, -1, 1, 1, 1, 0, 1) 
-5(0,0,0,1,0,1,1) -5(0,0,0,1,1,0,1) 
-25(0,0,1,1,0,0,1) +5(0,0,1,1,0,1,0) 
+25(0, 0, 1, 1, 1, -1, 1) - 5(0, 0, 1, 1, 1, 0, 0) 
+5(0, 1, -1, 1,0, 1, 1) + 25(0, 1, 0, 1, -1, 1, 1) 
-25(0, 1,0, 1, 0, 0, 1) - 5(0, 1,0, 1,0, 1,0) 

+5(0,1,0,1,1,0,0)] (B.l) 

CT), 2-loop = M. ( «\ r 64m 6 y2jB(1) j x x ^ Q) 1} + 48m 4 y2i?(0) ^ ^ ^ ^ Q ^ 

-32m 4 y 2 B(l, 0, 1, 1, 1, 0, 1) - 48m 4 y 2 5(l, 1, 1, 0, 1, 0, 1) 
+32mV5(l, 1, 1, 1, 0, 0, 1) + 16?nV5(l, 1, 1, 1, 1, 0, 0) 
-32m 4 y5(l, 0, 1, 1, 1, 0, 1) - 32m 4 y5(l, 1, 0, 1, 1, 0, 1) 
-12m 2 y5(0, 0, 1, 1, 1, 0, 1) - 12m 2 yB(0, 1, 0, 1, 1, 0, 1) 
-12m 2 y5(0, 1, 1, 1, 0, 0, 1) - 12m 2 y5(0, 1, 1, 1, 1, -1, 1) 
+24m 2 y5(0, 1, 1, 1, 1, 0, 0) + 16m 2 y5(l, -1, 1, 1, 1, 0, 1) 
+ 16m 2 y5(l, 0, 0, 1, 1, 0, 1) + 24m 2 y5(l, 0, 1, 0, 1, 0, 1) 
-16m 2 y5(l, 0, 1, 1,0, 0, 1) + 4m 2 y5(l, 0, 1, 1, 1, -1, 1) 
-12m 2 yB(l, 0, 1, 1, 1, 0, 0) + 24m 2 y5(l, 1,0, 0, 1, 0, 1) 
-12m 2 y5(l,l, 0,1,0,0,1) - 12m 2 y5(l, 1, 0, 1, 1, 0, 0) 
-8m 2 y5(l, 1, 1, 1, 0, -1,1)+ 8m 2 yB(l, 1, 1, 1, 0, 0, 0) 
+8m 2 y5(l, 1, 1, 1, 1, -1, 0) - 8m 2 y5(l, 1, 1, 1, 1, 0, -1) 
+4m 2 5(l, -1, 1, 1, 1, 0, 1) - 8m 2 5(l, 0, 0, 1, 1, 0, 1) 
+4m 2 5(l, 1, -1, 1, 1, 0, 1) + 35(0, 0, 1, 1, 0, 0, 1) 
-35(0, 0, 1, 1, 1, -1,1)- 35(0, 1, 0, 1, 0, 0, 1) + 35(0, 1, 0, 1, 1, -1, 1) 
-25(1,-2,1,1,1,0,1) +45(1,-1,0,1,1,0,1) -35(1,-1,1,0,1,0,1) 
+25(1,-1,1,1,0,0,1) -5(1,-1,1,1,1,-1,1) +25(1,-1,1,1,1,0,0) 
-25(1,0,-1,1,1,0,1) +65(1,0,0,0,1,0,1) - 35(1,0,0,1,0,0,1) 
+5(1,0,0,1,1,-1,1) -45(1,0,0,1,1,0,0) +25(1,0,1,1,0,-1,1) 
-25(1,0,1,1,0,0,0) -25(1,0,1,1,1,-2,1) +25(1,0,1,1,1,-1,0) 
-35(1,1,-1,0,1,0,1) +5(1,1,-1,1,0,0,1) +25(1,1,-1,1,1,0,0) 
-25(1, 1,0, 1, -1,0, 1) + 25(1, 1,0, 1,0, -1, 1) 

+25(1,1,0,1,0,0,0) -25(1,1,0,1,1,-1,0)] (B.2) 
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7 



(3), 2-loop 




-4e 



x [32m 6 y 2 B(2, 1, 1, 1, 0, 0, 1) + 16m 4 y 2 B(l, 1, 1, 1, 0, 0, 1) 
+ 16m 4 y 2 B(2, 1, 1, 0, 0, 0, 1) - 16?nV 5(2, 1, 1, 1,0, 0, 0) 
-16m 4 y5(2, 0, 1, 1, 0, 0, 1) - 16m 4 y5(2, 1, 0, 1, 0, 0, 1) 
+ 16m 2 y5(0, 1, 1, 1, 0, 0, 1) - 4m 2 yB(l, 0, 1, 1,0, 0, 1) 
-4m 2 yB(l, 1, 0, 1, 0, 0, 1) + 4m 2 y5(l, 1, 1, 1, -1, 0, 1) 
+4m 2 yB(l, 1, 1, 1, 0, -1,1)- 8m 2 yB(l, 1, 1, 1, 0, 0, 0) 
-8m 2 yB(2, 0, 1, 0, 0, 0, 1) + 8m 2 yB(2, 0, 1, 1, 0, 0, 0) 
-8m 2 yB(2, 1, 0, 0, 0, 0, 1) + 8m 2 yB(2, 1, 0, 1, 0, 0, 0) 
+2m 2 B(2, -1, 1, 1, 0, 0, 1) - 4m 2 5(2, 0, 0, 1, 0, 0, 1) 
+2m 2 B(2, 1,-1, 1,0, 0,1) - 5(1,0,1,1,-1,0,1) 
+5(1, 0, 1, 1, 0, -1, 1) + B(l, 1, 0, 1, -1,0, 1) 
-B(l, 1,0, 1,0, -1, 1) + B(2, -1, 1,0,0,0, 1) 
-5(2,-1,1,1,0,0,0) - 25(2,0,0,0,0,0,1) 
+25(2, 0, 0, 1, 0, 0, 0) + 5(2, 1, -1,0, 0, 0, 1) 

-5(2,1,-1,1,0,0,0)] (B.3) 



x [32mV5(l, 2, 1,0, 1,0, 1) + 16m 4 y 2 B(l, 1, 1,0, 1, 0, 1) 
-16mV 5(1, 2, 1, 0, 0, 0, 1) + 16?n 4 y 2 5(l, 2, 1, 0, 1, 0, 0) 
-16m 4 y5(l, 1, 1,0, 1,0, 1) - 16m 4 y5(l, 2, 0, 0, 1, 0, 1) 
-8m 2 yB(l, 0, 1, 0, 1, 0, 1) + 8m 2 yB(l, 1, 1, 0, 0, 0, 1) 
-8m 2 yB(l, 1, 1, 0, 1, -1, 1) + 8m 2 yB(l, 2, 1, 0, 0, -1, 1) 
-8m 2 yB(l, 2, 1, 0, 0, 0, 0) - 8m 2 yB(l, 2, 1, 0, 1, -1,0) 
+8m 2 yB(l, 2, 1, 0, 1,0, -1) + 2m 2 5(l, 0, 1, 0, 1, 0, 1) 
-4m 2 5(l, 1, 0, 0, 1, 0, 1) + 2m 2 5(l, 2, -1, 0, 1, 0, 1) 
+5(1, -1, 1,0, 1,0, 1) - 25(1,0,0,0, 1,0, 1) 
-5(1,0,1,0,0,0,1) +25(1,0,1,0,1,-1,1) 
-5(1,0,1,0,1,0,0) +5(1,1,-1,0,1,0,1) 
-25(1,1,0,0,1,-1,1) +25(1,1,0,0,1,0,0) 
-25(1,1,1,0,0,-1,1) +25(1,1,1,0,0,0,0) 
+25(1, 1, 1,0, 1, -2, 1) - 25(1, 1, 1,0, 1, -1,0) 
+5(1,2,-1,0,0,0,1) -5(1,2,-1,0,1,0,0) 
+25(1,2,0,0,-1,0,1) - 25(1,2,0,0,0,-1,1) 

-25(1, 2, 0, 0, 0, 0, 0) + 25(1, 2, 0, 0, 1, -1,0)] (B.4) 
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C The IBP identities 



Here we list the IBP identities 



d d k d d l d 



f \ ACM)* 



o 



(2vr) a! (2vr) d dp^ 

where p = k,l and q = k, I, q±. We present these identities in the form 

0(p,q)B(m, . . . ,n 7 ) = 



(C.l) 



(C.2) 



where the operators 0(p,q) corresponding to the insertion of d/dp^q^ into the loop integral can 
be rewritten in terms of the operators (6.10) as 



Q(k,q. 



2m 2 n 5 5 + + 2m 2 n 7 7 + - n 4 + n 5 - 2m 2 (2y - l)n 6 6 + + n 7 7 + l" 
-n 7 7 + 2" - n 5 5 + 4 _ - n 6 6+4 _ - n 7 7 + 4" + n 4 4+5" + n 6 6 + 5~ + n 7 7 + 5" 



0(k, q+) = 2m 2 n 6 6 + + 2m 2 n 7 7 + - n 4 + n 6 - 2m 2 (2y - l)n 5 5 + + n 7 7 + l" 

-?i 7 7 + 3~ - n 5 5 + 4~ - n 6 6+4~ - n 7 7+4~ + n 4 4+6~ + n 5 5+6- + n 7 7 + 6~ 

0(1, q-) = 2m 2 n 1 l + + 2m 2 n 3 3 + + 2m 2 n 7 7 + + n x - n 2 + 2m 2 (2y - l)n 3 3 + + n 2 2+l~ 
+n 3 3 + l~ + n 7 7 + l~ - nil+2~ - n 3 3+2~ - n 7 7+2~ - n 7 7+4" + n 7 7 + 5~ 

0(1, q + ) = 2m 2 n 1 l + + 2m 2 n 2 2+ + 2m 2 n 7 7 + + n x - n 3 + 2m 2 (2y - l)n 2 2+ + n 2 2+l _ 
+?i 3 3 + l~ + n 7 7 + l- - ml + 3~ - n 2 2+3 _ - n 7 7+3~ - n 7 7+4~ + n 7 7+5~ 

0(k, I) = 2m 2 n 5 5 + + 2m 2 n 6 6 + + 2m 2 n 7 7 + + n 7 - n 4 - n 4 4 + l" + n 7 7 + l~ - n 5 5 + 2" 
-?i 6 6 + 3~ - n 5 5 + 4" - n 6 6 + 4~ - n 7 7+4" + n 4 4+7" + n 5 5+7~ + n 6 6+7~ 

0(Z, fe) = -nil + 4~ + nil + 7~ - m + n 7 - n 2 2+l~ - n 3 3 + l~ - n 7 7+l~ 
+n 7 7 + 4" - n 2 2+5~ - n 3 3 + 6~ + n 2 2 + 7~ + n 3 3+7~ 

0(/, I) = d — 2n\ — n 2 — n 3 — n 7 — 2rr?n\\ Jr — 2m 2 n 7 7 + 
-n 2 2 + l" - n 3 3+l~ - n 7 7+l" + n 7 7+4~ 



Q(k, k) = d - 2n 4 - n 5 - n 6 - n 7 + n 7 7 + l - n 5 5 + 4 - n 6 6 + 4 - n 7 7 + 4 . 



(C.3) 



D Results for Master Integrals 

The results for MI are written in terms of the dimensionless quantities b(n\, . . . ,n 7 ) defined in 
(6.17). 
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D.l Two propagator topology 

This topology contains only one MI depicted in Fig. 14. It factorizes to one loop diagrams (i.e. a 
square of a tadpole) and hence can be computed straightforwardly. We get the result 



6 = 6(0,0,0,0,0, 1,1) = 




Figure 14: The two propagator MI. 



which has the following expansion 



r 2 ^ 1 ) (di) 
r 2 ( e + i)' { } 



&(- 2 ) = 1 

= 2 
&(°> = 3 
&« = 4 

= 5 . 



D.2 Three propagator topology, type a 

There is a single MI in this topology, see Fig. 15 . It is y independent thus can not be computed 
using the differential equations technique as formulated in this paper. Moreover, it does not 
factorize so the evaluation is more involved. However, the result can be found in the literature. 
The pioneering calculation was done in [38] as far as we know. It was further generalized in [39] 
and both results agree. The expansion in powers of e reads [38, 39] 



b = 6(1, 0,0, 0,0, 1,1) = 




Figure 15: The three propagator MI, type a. 



b (-2) 


3 

~ 2 


6 (-D 


17 

" T 


6 (0) 


59 

~ 1T 


6 (1) 


/ 65 

= 32 

V512 


h m 


= -64 f— 



+ 



7T 

24 



7C(3) + 1117 



16 



2048 



13n 2 
"96" 



1 



7T 2 l0g(2) 
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D.3 Three propagator topology, type b 

This single Master integral depicted in Fig. 16 does not depend on y and could be calculated 
directly using Feynman parametrization. 



b = 6(0,0, 1,1,0,0, 1) = 




Figure 16: The three propagator MI, type b. 



T(3 - 4e)r(l - e) 2 r(e)r(2e - 1) 
r(3-3e)r(2-2e)r(e + l) 2 



(D.2) 



Expanding the above formula gives 

6 (-2) 



6 (-D 
b (0) 
b (l) 
6 (2) 



24 ( 33 + 87r2 ) 



— (26 + 10vr 2 + 5^ (2) (1) + 8V> (2) (2) 



12 



— (660vr 2 + 256vr 4 + 15(-751 + 5O0 (2) (1) + 80i/> (2) (2) - 37O0 (2) (3)) 



D.4 Three propagator topology, type c 

This topology includes one MI (see Fig. 17) and corresponds to the product of one-loop integrals. 
The e expansion of the latter was studied in [40]. The convenient point for fixing the integration 



b = 6(0,0,0,0, 1,1,1) 




Figure 17: Three propagator MI, type c. 
constants was in this case x = 1 which is below the physical threshold. The solution is 
6(" 2 ) = 1 



6 (-i) 



(x + l)H(0,x) -3x + 3 
x — 1 



4(0) _ _3(x+l)ff(0,x) \ MgW + I f^ + h + 42 

x — 1 \x — 1 / 1 — x 6 \ x — 1 

= ^dfW-M) (^-42)(x + l)fl(0,x) N 

3(x-l) 6(x-l) V^- 1 / 

6(x + l)ff(-l,0,x) 3(x + l)ff(0,0,x) 4(x + l)ff(-l, -1, 0, x) 

x— 1 x — 1 x — 1 

+ ( 4 + A ff Q) Qj ^ + (x + l)g(0,0,0,x) + 4(x + 1)C(3) + 30(x - 1) + vr 2 (x + 1) 



x - 1 J x ' ' ' ' 1 - x 2(x - 1) 



54 



We checked this MI with [43] where it is computed using the same formalism to 0(e 1 ) including. 
Comparing the results gives an agreement. 

D.5 Three propagator topology, type d 

This MI depicted in Fig. 18 could be computed directly using loop integration (it factorizes to an 
off-shell massless bubble and a tadpole). The result reads 



6 = 6(0, 1,1,0,0,0, 1 



Figure 18: Three propagator MI, type d. 

4-H/)-T(i-6) 2 r(6-i)r( g ) 

r(2-2e)r(e + l) 2 

and has the following expansion 
6(" 2 ) = 1 

ftC" 1 ) = H(0,x) +2H(l,x) + 3 




(D.3) 



&(°) = SH(0, x) + 6H(l,x) + 2H(2, x) + H(0, 0, x) + 2H(1, 0, x) + 4H(1, l,x) + 7 

6 

+12if (1, 1, x) + 4ff(l, 2, x) + 2#(2, 0, x) + 4ff(2, 1, x) + H(0, 0, 0, x) + 2H(1, 0, 0, x) 



^ 2 



+4H"(1, 1,0, x) + 8#(1, 1, 1, x) - 2C(3) - y + 15 



b {2) = (-2C(3) + 15 - y) H(0, x) + (-4C(3) + 30 - vr 2 ) #(1, x) + (l4 - #(2, *) 



vr 2 \ „, , / vr 2 



+6ff(3, x) + 2#(4, x) + ( 7 - -) fl-(0, 0, x) + Il4 - —j H(l, 0, x) 

+ ^28 - ^\ H(l, 1, x) + 12H(1, 2, x) + 4ff(l, 3, x) + 6if(2, 0, x) + 12iT(2, 1, x) 

+4F(2, 2, x) + 2H(3, 0, x) + 4ff(3, 1, x) + 3#(0, 0, 0, x) + 6#(1, 0, 0, x) 

+12if (1, 1, 0, x) + 24ff (1, 1, 1, x) + 8iT(l, 1, 2, x) + 4ff(l, 2, 0, x) + 8H(1, 2, 1, x) 

+2H{2, 0, 0, x) + 4H(2, 1, 0, x) + 8H(2, 1, 1, x) + £T(0, 0, 0, 0, x) + 2H(1, 0, 0, 0, x) 

vr 4 7tt 2 

+4ff(l, 1,0, 0, x) + 8tf (1, 1, 1,0, x) + 16#(1, 1, 1, 1, x) - 6C(3) - — - — + 31. 

4U b 

Also this MI was checked with [42] and a complete match was reached. 
D.6 Three propagator topology, type e 

There are two nested Mis (see Fig. 19). The suitable point for fixing the integration constants 
was chosen = 1. Then the solution reads 

(_i) _ xH(0,x) 
1 ~ x 2 -l 
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Figure 19: The two nested Mis with three propagators; type e. 





QxH(-l,0,x) 2xH(l,0,x) (5x - 3)xH(0, 0, x) n 2 x 

x 2 — 1 x 2 — 1 (x — l) 2 (x + 1) 6 — 6x 2 

2x#(0,0,x) 
(x - l) 2 

tt 2 xH(-1,x) tt 2 xH(1,x) 36xff(-l,-l,0,x) 24xH(-l, 0, 0, x) 

x 2 — 1 3 — 3x 2 x 2 — 1 x 2 — 1 

_12xiT(-l,l,0,x) 12xif(l,-l,0,x) 4x#(l,l,0,x) 
x 2 — 1 x 2 — 1 x 2 — 1 

tt 2 (3- 5x)xff(0,x) 6(3- 5x)xff(-2,0,x) 2(5x - 3)xff(2,0,x) 
+ 6(x-l) 2 (x + l) + (x- l) 2 (x + l) + (x-l) 2 (x + l) 
(13a - 7)xff(0,0,0,x) 2(3x - 5)xff (1, 0, 0, x) 2(4 - 7x)xC(3) 
+ (x- l) 2 (x + l) + (x- l) 2 (x + l) + (x-l) 2 (x + l) 
vr 2 xi7(0,x) 12xtf(-2,0,x) 4x#(2,0,x) 6xH(0, 0, 0, x) 
3(x-l) 2 + (x- l) 2 (x - l) 2 (x - l) 2 

4x#(l,0,0,x) 6xC(3) 

+ (x - I) 2 + (x- l) 2 

tt 2 (5x - 3) J H"(-2,x)x vr 2 (3 - 5x)H(2,x)x 6(7 - 13x)if (-3, 0, x)x 
(x-l) 2 (x + l) + 3(x - l) 2 (x + 1) + (x - l) 2 (x + 1) 
6ir 2 H(-l,-l,x)x 47r 2 jj'(-l,0,x)x 2vr 2 ff(-l,l,x)x 

X 2 — 1 X 2 — 1 X 2 — 1 

vr 2 (7- 13x) J H"(0,0,x)x 2vr 2 #(l,-l,x)x vr 2 (5 - 3x)H(l, 0, x)x 
+ 6(x - l) 2 (x + 1) h x 2 - 1 + 3(x - l) 2 (x + 1) 

27r 2 ff(l,l,x)x 2(13x - 7)ff(3,0,x)x 36(5x - 3)H{-2, -1, 0, x)x 
+ 3-3x 2 + (x- l) 2 (x + l) + (x-l) 2 (x + l) 

24(3 - 5x)ff(-2,0,0,x)x 12(3 - 5x)H (-2, 1, 0, x)x 144ff(-l, -2, 0, x)x 

+ (X-I) 2 (x+1) + (X-I) 2 (x + 1) + X^l 

48ff(-l,2,0,x)x 12(5 - 3x)ff(l,-2,0,x)x 4(3x - 5)^(1, 2, 0, x)x 

x 2 - 1 + (x - l) 2 (x + 1) (x - l) 2 (x + 1) 

12(3 - 5x)ff(2,-l,0,x)x 2(27x - 17)ff(2, 0, 0, x)x 4(5x - 3)H(2, 1, 0, x)x 

+ (X-1) 2 (X+1) + (X- l) 2 (x + l) + (X- l) 2 (x + l) 

216ff(-l,-l,-l,0,x)x 144ff(-l,-l,0,0,x)x 72H(-1, -1, 1, 0, x)x 
X 2 — 1 x 2 — 1 x 2 — 1 



5G 



60JT(— 1, 0, 0, 0, sc)a; 72iJ(-l, 1, -1, 0, x)x 48iJ(-l, 1, 0, 0, x)x 

x 2 — 1 x 2 — 1 x 2 — 1 

24ff(-l,l,l,0,x)x (29x - 15)g(0,0,0,0,g)g 72ff(l, -1, -1, 0, x)x 

x 2 — 1 (x — l) 2 (x + 1) x 2 — 1 

48ff(l,-l,0,0,x)x _ 24ff(l,-l,l,0,x)x 2(7x - 13)ff(l, 0, 0, 0, x)x 

x 2 — 1 x 2 — 1 (x — l) 2 (x + 1) 

24ff(l,l,-l,0,x)x 4(5x -3)ff(l, 1,0,0, x)x 8ff(l, 1, 1, 0, x)x 

x 2 — 1 (x — l) 2 (x + 1) x 2 — 1 

tt 4 (35 - 61x)x 66i?(-l,x)C(3)x (18 - 34x)#(0, x)C(3)x 

"^TTTTTTT TToT ! TT "I - o ^ 



360(x - l) 2 (x + 1) x 2 -l (x-l) 2 (x + l) 

4(7 - 4x)H(l, x)C(3)x 
(x-l) 2 (x + l) 

{2 ) _ 16xC(3)iJ(0,x) 12xC(3)#(l,x) 2vr 2 xi?(-2,x) 2tt 2 xH(2,: 



+ 



(x-1) 2 (x-1) 2 (x-1) 2 3(x-l) 2 

36xiJ~(-3, 0, x) vr 2 xi?(0,0,x) 2vr 2 x#(l,0,x) 12x#(3,0,x) 



(x-1) 2 (x-1) 2 3(x-l) 2 (x-1) 2 

72xff(-2,-l,0,x) 48x^-2,0, 0,x) 24xff (-2, 1, 0, x) 

(x-1) 2 + (x-1) 2 + (x - l) 2 
24xF(l,-2,0,x) 8xiJ(l,2,0,x) 24x#(2, -1, 0, x) 20xH(2, 0, 0, x) 



(x-1) 2 (x-1) 2 (x-1) 2 (x-1) 2 

8xiT(2,l,0,x) UxH(0, 0, 0, 0, x) 12xff(l, 0, 0, 0, x) 8xff(l, 1, 0, 0, x) 13vr 4 x 

H~ 7 TTo 7 TTo ~T~ 



(x-1) 2 (x-1) 2 (x-1) 2 (x-1) 2 180(x-l) 2 ' 

This pair of Mis belongs to the general class J011, the e expansion of which was studied system- 
atically in [40, 41]. In the framework of differential equations approach they were calculated in 
different basis across the literature [43, 44, 45]. We use the same basis as in [44] thus we verified 
their results to 0(e 2 ) including. 

D.7 Four propagator topology, type a 

This topology includes the Master integral depicted in Fig. 20. It is a product of a massive (see 



6 = 6(0,1,1,0, 1,1,0) 



Figure 20: The four propagator MI, type a. 

[40]) and a massless bubble. The former one can be easily extracted from T3c and the latter one 
can be calculated directly. So we used this trick instead of solving the corresponding differential 
equation. The solution for current orders reads 




= 1 

2H(0,x) | 2 
x — 1 

4 \ 8H(0,x) 



x — 1 



b<® = 8iJ(l, x) + f — + 2 J H(-l, 0, x) + 4i?(l, 1, x) 



x - 1 
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4#(2, x) 2(x + 2)H(0, 0, x) 4#(1, 0, x) tt 2 ^ 



x- 1 



x — 1 



x-1 3(x-l) 



5 (D 



vr 2 (x + 1)#(-1, x) (vr 2 (x + 2) - 72) #(0, x) 



+ 



+ 



+ 



3(x - 1) 
2vr 2 



+ 



3(x - 1) 



+ 24 #(l,x) 



x — 1 



+ 4 ff(-2,0,i) + 



3(x - 1) 
16#(2,x) 4(x + 2)#(3,x) 
x — 1 x — 1 

(x + l)#(-l,0,x) 



x — 1 
+16H(l,l,x] 



+ 4 #(-l,2,x) 



x — 1 

(x + 2)#(0,0,x) 16#(l,0,x) 



x — 1 x — 1 

8#(l,2,x) 4(x + 2)H (2, 0, x) 8#(2,l,x) 



x — 1 



x- 1 



x — 1 



4(x + l)#(-l,-l,0,x) 6(x + l)#(-l,0,0,x 



x - 1 



+ 



x - 1 
14 



+ 



x- 1 



+ 4 #(1,-1,0,*) 



x- 1 
4(x + 2)#(l,0,0,x) 
x — 1 



8#(l,l,0,x) 4 (3(8x + £g) - 8) + tt 2 ) 
— ^ + 8#(1, 1, 1, x) + _ . 



Making a check with [42] up to 0(e 1 ) lead to full agreement. 



D.8 Four propagator topology, type b 



6 = 6(0, 1,0, 1,0, 1,1) 




Figure 21: The four propagator MI, type b. 

The solution for the one MI depicted in Fig. 21 is (the integration constants were fixed at 
x = 1) 

1 

2 
5 

2 



6 (-2) 
6 (-D 



( 0) _ 2vr 2 x#(0,x) _ 4x#(0,0,0,x) _ 19 
3-3x 2 x 2 -l ^lU,U,xj+ 2 



6 (D 



(24xC(3) + ^ 2 ((x - 4)x - 1)) #(0, x) 24x#(-3,0,x) 4tt 2 x#(-1, 0, x) 
6 (x 2 — 1) x 2 — 1 3 — 3x 2 

(-5x 2 + 2^ 2 x + 5) #(0,0, x) 4^ 2 x#(l,0,x) 8x#(3,0,x) 



x 2 - 1 3 (x 2 - 1) x 2 - 1 

(3-x(3x + 4))#(0,0,0,x) 8x#(2,0,0,x) 8x#(-l, 0, 0, 0, x) 



+- + 



x 2 — 1 x 2 — 1 x 2 — 1 
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4x#(0, 0,0,0, a) 8xtf(l,0, 0,0, x) 

H n ; 1- 6-H(— 2,0,xj 



x 2 — 1 x 2 — 1 



2#(2,0,x) 



+2H (1, 0, 0, x) + 



45 (x 2 - 1) (6C(3) + 65) + 26vr 4 x 



90(x 2 -l) 

A total correspondence with [43] to 0(e 1 ) was found. 
D.9 Four propagator topology, type c 

The third four propagator topology contains one Master integral depicted in Fig. 22 which has 



6 = 6(1,1,1,0, 0,0,1) 




Figure 22: The four propagator MI, type c. 

the following solution (integration constants fixed at x = —1) 
(_!) _ 2xH(2,x) , xH{0,0,x) , 2vr 2 x 



+ - + 



x 2 — 1 x 2 — 1 3 (x 2 — 1) 

(0 ) _ 4vr 2 x#(-l,x) t it 2 xH(0,x) t 4ir 2 xH(l,x) i 2xH(2,x) t 2xH(3,x) t 4xH(-l,2,x) 



+ — = + » , o x + o ' + o ' + 



3 (x 2 — 1) 6 — 6x 2 3 (x 2 — 1) x 2 — 1 x 2 — 1 x 2 — 1 

xH(0,0,x) 4xH(l,2,x) 2xH(2,0,x) 4xH(2,l,x) 2xH(-l, 0, 0, x) 

I o I Q I q j I f-i _ j I 



c 2 — 1 X 2 — 1 x 2 — 1 x 2 — 1 x 2 — 1 

xff(0,0,0,x) 2x^(1,0, 0,x) x (15C(3) + 2tt 2 ) 
x 2 — 1 x 2 — 1 3 (x 2 — 1) 

2x (15C(3) + 2vr 2 ) H(-l, x) 2x (l5((3) + 2vr 2 ) H(l , x) x (12£(3) + 7T 2 ) .ff(0, 



x 



3(x 2 -l) 3(x 2 -l) 6-6x 2 

(vr 2 -6) xH(2,x) 2xH(3,x) 2xH(4,x) 87r 2 xff(-l, -1, x) tv 2 xH(-1, 0, : 

I q _, I q _, I <-» / rj \ I 



3(x 2 -l) x 2 -l x 2 -l 3(x 2 -l) 3-3x 2 

8tt 2 x j H"(-1,1,x) AxH(-l,2,x) 4xx7(-l,3,x) (vr 2 - 6) xiJ(0, 0, x) 

3 (x 2 — 1) x 2 — 1 x 2 — 1 6 (x 2 — 1) 

87r 2 x#(l,-l,x) tt 2 xH(1,0,x) 8ir 2 xH(l,l,x) 4xH(l,2,x) 4xH(l, 3, x) 

I ~ / o -■ \ I . . o I r-« / n \ I o I 



3(x 2 -l) 3-3x 2 3(x 2 -l) x 2 -l x 2 -l 

2xff(2,0,x) 4xH(2,l,x) 4xH(2,2,x) 2xH(3,0,x) 4xff(3,l,x) 

x 2 — 1 x 2 — 1 x 2 — 1 x 2 — 1 x 2 — 1 

8xff(-l,-l,2,x) 2xff(-l,0,0,x) 8xiJ(-l,l,2,x) 4xiT(-l, 2, 0, x) 

I o I o .. I o I 



x 2 — 1 x 2 — 1 x 2 — 1 x 2 — 1 

8xff(-l,2, l,x) xff(0,0,0,x) 8xff(l, -1, 2, x) 2xff(l,0,0,x) 8xiT(l, 1, 2, x) 

I rt I r> I rt 1 rj I 



x 2 — 1 x 2 — 1 x 2 — 1 x 2 — 1 x 2 — 1 

4x#(l,2,0,x) 8xff(l,2,l,x) 2xF(2,0,0,x) 4xH(2, 1, 0, x) 8xF(2,l,l,x) 



x 2 — 1 x 2 — 1 x 2 — 1 x 2 — 1 x 2 — 1 

4xff(-l,-l,0,0,x) 2xff(-l,0,0,0,x) 4xff(-l,l,0,0,x) xH (0, 0, 0, 0, x) 

I rj -, I rj -i I o -t I 



c 2 — 1 x 2 — 1 x 2 — 1 x 2 — 1 

| 4xff(l,-l,0,0,x) 2xff(l,0,0,0,x) 4xff(l,l,0,0,x) + x (1800C(3) + 240^ 2 + 59tt 4 ) 
x 2 — 1 x 2 — 1 x 2 — 1 360 (x 2 — 1) 
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We verified the result of [42] to 0(e 1 ). 

D.10 Four propagator topology, type d 

Our choice for the two Mater integrals in this topology is depicted in Fig. 23. Here we tried two 



61 = 6(1, 1,0,0,0,1,1) = 



& 2 = 6(2,1,0,0,0,1,1) = 




Figure 23: The basis of MI for the four propagator topology, type d. 

points x = ±1 for fixing the integration constants and both gave the same result. The expansion 
of the solution in powers of e is 



b (-2) 



6<-» 
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2 
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2 
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2 



( _D = (a + l)ff(0,a) 
2 2-2x 



b 



(0) 2ir 2 xH(0,x) 4aff(0,0,0,a) 1 2 

1 = ^3^2 ^(0,0,x) + -(57-2 7 r) 
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| 6(x + l)ff(-l,0,x) 



6(x - 1) 
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x — 1 



x — 1 



x — 1 / 
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These Master integrals appear in [43] in the same basis. Making a check leads to complete agree- 
ment up to 0(e 1 ). 



D.ll Four propagator topology, type e 

This topology includes three coupled Master integrals 6( — 1, 1, 1, 1, 0, 0, 1), 6(0,1,1,1,0,-1,1), 
6(0,1,1,1,0,0,1). We make a transformation to a more suitable basis depicted in Fig. 24. In 
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62 = 6(0, 1,1,2,0,0,1) 




6 3 = 6(0, 1,1, 1,0, 0,1) = 




Figure 24: New basis of Mis for the four propagator topology, type e. 



this basis the system decouples and the solution may be found as in previous cases (using the 
point x = —1 to fix the unknown integration constants) 
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(0) (2x 3 -2x + 4)g(2,x) (x 3 - x + 2) H (0, 0, x) (7x 7 \ 
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The last two Mis 62 and ^3 can be found in [42] in the same basis while for the first one we use a 
different basis. The former two were checked to 0(e 1 ) and the latter one was left uncontrolled. 

D.12 Five propagator topology, type a 

The current topology is formed by the MI depicted in Fig. 25. It corresponds to F10101 of [40, 41] 
where the e expansion was given. This MI has the following expansion using x = 1 for fixing the 



6 = 6(0, 1,1,0, 1,1,1) 




Figure 25: The five propagator MI, type a. 
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integration constants 
6(°) = 
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Verifying the result with [42] to 0(e 1 ) gave full agreement. 



D.13 Five propagator topology, type b 

We have chosen the two Master integrals as depicted in Fig. 26 (the AIR original basis was 
6(0,1,1,1,0,1,1) and 6(0,1,1,1,-1,1,1)). This choice has been done in order to decouple the 
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6l = b(-l,l, 1,1, 0,1,1) = X 



6 2 = 6(0, 1,1, 1,-1, 1,1) 




Figure 26: Suitable basis of Mis for the five propagator topology, type 6. 

system of corresponding differential equations. The appropriate point for fixing the integration 
constants was a = — 1. Then it has the solution 



1 

2 
1 

2 
5 

2 

6^ 1} = H(0,x) + 2H(l,x) + ^ 



b[~ 2) 



,(-2) 
°2 



(0) _ ^^^^^ 

bl - a2-l + 



3x 2 ff(0,0,0,x) (vr 2 (3a 2 - l) + 12a(a + 1)) H (0, a) 2H(3, 



6 (a 2 — 1) a 2 — 1 



64 



[x (-2x 2 + x - 2) - l) H(0, 0, x) + (x 2 + l) H(2, 0, x) 



(x-l) 2 (x + l) 1-x 2 
- ''^y + 2IT<1, ,) + HS^ + H(l. 0, x) + 2* (1, l,x) 

o'o. x) + 12(57 - 12C(3)) + 3 t - 11 + 24C(3) - 57 

6 (x^ — 1) 

(3(x - l)(4(x - 1)C(3) + 5x + 5) + vr 2 (x(3x - 1) - 1)) H(0, x) 

3(x 2 -l) 

2(x - 2)H(3,x) 2((x - l)x + l)#(2,0,x) 4((x - l)x + 1)H(2, 1, x) 



x 2 — 1 x 2 — 1 x 2 — 1 

3x(2x-l)ff(0,0,0,x) ^ 2(x-l)C(3) 
+ x^l + I x + 1 + 10 I H{hX 



+ (^TT- 2 )^^ + (xTT- 2 )^^ 



(- (24 + 7T 2 ) x + vr 2 + 6) F(0, 0, x) (vr 2 (x - 1) + 9(x + 1)) £T(1, 0, x) 



6(x + l) 3(x + l) 

-6H(1, 1, x) + f — ^ - 2^ (1, 3, x) + f - H(3, 0, 

2 J #(3,l,x) + 4ff(l,0,0,x) 



x + 1 



+ (-JTi " 2 ) ff 2 ' °> x ) + (t+T " 4 ) H (1 ' 2 ' X ' X) 

+ " 2 ) F ( 2 ' °' °' x ) + ( 3 " ^l) H ^ °' °' °' x ) 

-72 (2x 2 + x - 4) C(3) + 342 (x 2 - l) + vr 4 (x - l) 2 - 12vr 2 (7x - 3)(x - 1) 



36 (x 2 - 1) 



ff(0,0,0,0,x)x 2 2 2 (l2(x-l)(x + 6) + vr 2 (9x 2 -7))ff(2,x) 

2 — i h q 71 " -^l -2 ^) H ~eTT~2 — 7^ 

x z — 1 3 6 [x z — 1) 

2(x(x + 3)-5)ff(3,x) 2 (x 2 - 4) g(4, x) _ 6 (3x 2 + l) H{-3, 0, x) 

x 2 — 1 x 2 — 1 x 2 — 1 

12((x-l)x + l)if(-2,0,x) 

H ? -Tg h 2H(-2,2,x) 

[x — l) z 

^ (tt 2 (x 2 + 1) - 18(x + l) 2 ) H(-1,0, x) 



3 (x 2 - 1) 



+ 



/ 4 \ , . 2(x 2 + l)#(-l,3,x) 

— - - 2 -l,2,i - [ \ 

\x + 1 J x z — 1 

(vr 2 (x - 1) (lOx 2 + 3) + 6(x(x(2x - 9) + 14) + 5)) H(0, 0, 



6(x-l) 2 (x + l) 

(2 - 2x(2(x - l)x + 5))H(2, 0, x) 



(3(x + l)(9x - 5) - vr 2 (x 2 + 3)) H(l, 0, x) 

+— ' 2 \, " + 14H(1, 1, x 

3 [x z — 1) 



+ [l + — l )H{l,2, X ) + 2H( 1 , 3 , x)+ (X _ 1)2(I + 1) 

4(2x -3)ff(2,l,x) 2(x 2 + 5)ff(2,2,x) (7x 2 - 3) F(3, 0, x) 



X 2 — 1 X 2 — 1 X 2 — 1 



65 



2(x 2 -5)H(3, l,x) , , / 2 

+— 2 1 - + g(-2,0,0,a;)+ — — -1 g(-l,0,0,a;) 

x z — 1 \x + 1 

2 (x 2 + l) H(-l, 2, 0, x) 4 (x 2 + l) H(-l, 2, 1, 



X 



X 2 — 1 X 2 — 1 

(x 3 - lOx - 3) i/(0,0,0,x) 
(x — l) z (x + 1) 

+ + 9 + ^ + (a / 1)2 ) °> °> x ) + ^t 1 ' X > °> x ) 

, (5x 2 + 9)i?(2,0,0,x) 

+12iJ(l, 1, 1, x) + 4ff(l, 2, 0, x) + ^ ^ — V ' 

1 — x z 

6 (x 2 + l) H(2, 1, 0, x) 12 (x 2 + l) H{2, 1, 1, x) 3 (x 2 + l) H(-l, 0, 0, 0, : 



x 2 — 1 x 2 — 1 x 2 — 1 

(5x 2 - 13) iT(l, 0, 0, 0, x) iT(l, x) (vr 2 (7 - 9x) + 12(x + 1)(7 - 6C(3))) 



+ 



+ 



x 2 - 1 6(x + 1) 

H(-l,x) (6 (x 2 + 1) C(3) - 2^ 2 (x - l) 2 ) 
3 (x 2 - 1) 

H(0, x) (84x (x 2 - 1) + tt 2 (x + l)(x(10x - 13) + 5) - 6(x - 1) (4x 2 - 3) C(3)) 

6(x - l) 2 (x + 1) 

(-2925(x + l)(x - l) 2 + vr 4 (27x 2 - 17) (x - 1) 



90(x - l) 2 (x + 1) 
+ 15vr 2 (x(17x - 18) + 5)(x - 1) + 90(x(x(9x - 22) - 3) + 4)C(3)) 



2tt 2 (x - l)H(-3, x) (7vr 2 (x - l) 2 + 12(x - 2)(4x - 5)) H(3, 



3(3 + 1) 



6(x 2 - 1) 



+ 
+ 
+ 



1 



x + 1 

12 
x + 1 

4 



+ 



H(A,x 



x + 1 
+4H(-2,2,x) -f 



x - 1 
6 J H(-4,0,x) - 

2] H(-3,2,x) + 
4 



16 



8 H(5,x) 



x + 1 

12(x(3x - 2) + l)g(-3,0,a) 
x 2 — 1 

(vr 2 (x - 1) + 54(x + 1)) H(-2, 0, : 



x + 1 

4((x-l)x + l)#(-l,3,x; 
x 2 — 1 



2 ff(-2,3,x) + 



3(x + l) 
20 



10) #(-1,2, 
x + 1 / 



+ 



X + 1 J 



4 )H(-l,4,x) 



^ (vr 2 (x - 1) - 30(x - 3)) H(l, 2, x) + 8z£T(l, 3, x) 



+ 



20 

x + 1 



3(x + l) 
10 J #(1,4, a;) + 



17 + 2vr 2 
x + 1 



x + 1 



TT 2 - 10 + 



1 — X 



#(2,0,*) 



4(2(x - 4)x + 7)H(2, 1, x) 4((x - 3)x + 5)H(2, 2, x) 



x 2 -l 



x 2 - 1 



+ 



4 2] H(1, 3, x) + &!& x ~^>~ i ) H ( s < °> x ) 



X + 1 J 



| 8(x 2 -2)ff(3,l,x) 
x 2 — 1 



20 
x + 1 



x 2 — 1 
10] #(3,2, a?) + 



x + 1 J 



3) #(4,0,: 



66 



+ ( — — -r - 10 ) H(4, 1, x) + ( — - 1 ) F(-3, 0, 0, x) + 2#(-2, 0, 0, x) 
\x + 1 J \x + 1 / 

+ (^T " 2 ) F( ~ 2 ' 2 ' °' x) + (xTT " 4 ) H{ ~ 2, 2 ' x) 

(15 + ^ 2 ) (x - !)#(-!, 0, 0, x) + 2ir 2 (x - 1)H(-1, 1, 0, x) 



3(x + l) 3(x + l) 

4((x- l)x + !)#(-!, 2, 0, x) 



8((* -1)* + !)*(-!, 2, + ( 4 N ^ 



x 2 — 1 \X + 1 

(vr 2 (x - l) 2 + 26x 2 - 50x + 2) #(0, 0, 0, x) 



+ (^_-4)tf(-l,3,l,x 



2 (x 2 - 1) 



48 \ , . . 2tt 2 (x- l)i?(l,-l,0,x) 

+ - 24 H{1, -3, 0, x) - 24#(1, -2, 0, x) + 



x + 1 "7"^' 1 3(x + l) 

+ (l0 - 2 3 ( ^ + ~ 1) 1) ) #(1. 1. 0> + 20ff(l, 1, 1, x) + ( _ 4) 1, 3, x) 

( 8 \ , x 8(x- l)-ff(l,2, l,x) / 24 \ , 

+ 12 - — — H(l, 2, 0, x + > \' ' ' > + — — - 12 if 1, 2, 2, x 

\ x + iy x + i V x + i / 

( 4 \ . N / 24 \ , x 12(x - 1)^(2, -2,0, x) 

+ ( 2 - ) tf(l,3,0,x) + ( — ^ - 12 ) #(1,3,1,*) + 

+ 

x z — I 

24((x-l)x + l)ff(2,l,l,x) / .\ ff , 9n / 18 Q \ mqnn 

+ {J^i " 6 ) F(3 ' lj °' x) + (xTT " 12 ) F(3 ' 15 lj x) + ( 3 ~ xTi) f( ~ 2 ' °' °' °' x) 

6((x - l)x + !)#(-!, 0,0,0,x) 

x 2 — 1 V x + 1 



x + 1 / \x + 1 y x + 1 

(-6(x- l)x - 14)ff(2,0,0,x) _ 12((x- l)x + 1)H(2, 1, 0, x) 

x 2 — 1 x 2 — 1 

24((x-l)x + l)iT(2,l,l,x) , / 8 A / 18 



ti^py + » _ 4 ! 2 } 

x z — 1 \x + 1 / 

/ 16 \ , N / 8 \ , , x(2x - l)iT(0, 0,0,0, x) 

+ (^TT -4 ) ff ( 1 '- 1 > 2 .°' x ) + (jrp[~ 8 ) -1.2.1,1) 

+ 4(x(x + S )-g)H(l,0,0,0,x) + _ 4 j , 2, 0, i) + (j^j ~ 8j B(l, 1, 2, 1, x; 

+ (jr^J ~ 1S ) 2 .C °. + (j^TT " l 2 ) 2 . 1.0. + ( J^y - 24 j H(l, 2, 1, 1, : 
+ C » - ^ H (2,0,0,0,x) + «(*-l)ff(-U.O.O.O.*) + «(» -W, -1.0, 0.0,x) 



x + 1 y x + 1 x + 1 

+ (x-l)H(lAO,0,0,z) + (4 _ \ ^ A Qi 



+ 



X+1 \x + 1 

2fT(-l, 0, x) (l2C(3)(x - l) 2 + vr 2 ((x - l)x + 1)) 
3 (x 2 - 1) 



67 



#(1,0, x) (21C(3)(x - l) 2 - 4vr 2 ((x - 2)x + 2) + 45 (x 2 - l)) 
+ 3(x 2 -l) 

+H(-2, x) ^ — x +—)+ H(l, 1, x) ^ x + 1 + 30 j 

H(l,x) (3vr 2 (37 - 43x) - vr 4 (x - 1) + 684(x + 1) - 72(7x + 5)C(3)) 

18(x + 1) 

H(-l, x) (tt 2 (-60 + vr 2 ) (x - l) 2 + 72((x - l)x + 1)C(3)) 
+ 18 (x 2 - 1) 

H(0,x) (l7vr 4 (x - l) 2 + 1710 (x 2 - l) + 15vr 2 (x(2x + 13) - 5) - 90(x(16x - 17) + 2)C(3)) 
+ 90 (x 2 - 1) 

H(0, 0, x) (vr 2 (9(l - 2x)x - 4) + 6(x - l)(3C(3)(x - 1) - 8x + 5)) 
+ 6(x 2 - 1) 

H(2, x) (vr 2 (x(9x - 1) - 7) + 6(x - l)(6C(3)(s - 1) + x + 14)) 
+ 3(x 2 - 1) 

+ 7U7T^ TT (t 4 ((20 - 59x)x + 29) + 15vr 2 (x - l)(9((3)(x - 1) - 46x + 6) 

9U (x z — 1) 

+ 45 (-6C(5)(x - l) 2 + 65x 2 + (-44x 2 + 78x - 30) C(3) - 65)) 
| 4(x-l)ff(-l,l,x)C(3) | 4(x-l)ff(l,-l,x)C(3) 
x+1 x+1 
This couple of Mis can be found in [46]. We performed a transformation to the original basis 
and made a check only for 6(0, 1, 1, 1, 0, 1, 1) which exactly appears also in [46]. The result is in 
correspondence up to O(e ). However, this verifies both of our Mis since the transformation mixes 
them to 6(0, 1, 1, 1, 0, 1, 1). The above paper does not offer higher orders so we add them to the 
list of existing Mis. 
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